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ABSTRACT 


The  concepts  and  relations  which  form  the  theoretical 
foundation  of  the  statistical  energy  approach  to  vibration 
analysis  are  delineated.  The  utility  of  this  approach  for 
dealing  with  complex  systems  is  discussed,  and  its  range  of 
applicability  is  indicated. 

The  important  properties  of  the  modes  of  vibrating 
systems  are  reviewed  and  used  to  exhibit  relations  between 
modal  responses  and  average  responses  of  total  systems.  It 
is  demonstrated  that  under  some  conditions,  which  are  often 
approximated  in  actual  systems,  the  average  rate  of  flow  of 
energy  from  one  mode  to  another  is  proportional  to  the  dif¬ 
ference  in  the  modal  energies.  It  is  also  shown  that  under 
some  conditions  the  average  rate  of  flow  of  energy  between 
two  sets  of  modes  (representing  groups  of  modes  of  two 
coupled  systems  in  a  given  frequency  band)  is  proportional 
to  the  difference  in  the  set-average  modal  energies.  Ex¬ 
amples  are  presented  which  indicate  how  these  relations 
permit  one  to  obtain  approximate  solutions  to  complex  prob¬ 
lems  simply,  on  the  basis  of  energy  conservation  considera¬ 
tions. 

Available  extensions  of  the  major  concepts  developed 
in  detail  are  mentioned,  and  references  to  the  current 
literature  are  given. 
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INTRODUCTION 


Recent  years  have  seen  a  continuation  and  acceleration 
of  the  trend  toward  higher  performance  vehicles,  toward 
increased  propulsion  system  power,  and  toward  more  highly 
sophisticated  equipment  and  instrumentation.  This  trend 
has  brought  with  it  a  greater  incidence  of  problems  asso¬ 
ciated  with  high  frequency  vibrations  and  has  caused  aero¬ 
space  engineers  to  concern  themselves  increasingly  with 
vibrations  at  frequencies  considerably  above  the  fundamental 
structural  resonances. 

Classically,  vibration  engineers  have  focused  their 
attention  on  low-frequency  oscillations,  since  the  lowest 
few  vibration  modes  are  generally  the  ones  which  are  asso¬ 
ciated  with  the  greatest  deflections,  the  highest  stresses, 
and  gross  structural  failures.  However,  the  rather  complete 
arsenal  of  analysis  techniques  which  have  been  developed  for 
dealing  with  low-frequency  vibration  problems  contains  none 
which  can  deal  simply  and  effectively  with  most  high-frequency 
problems,  such  as  those  of  importance  in  relation  to  sonically 
induced  fatigue,  instrumentation  performance,  or  noise  trans¬ 
mission. 

Although  the  classical  methods  are  valid  in  principle 
at  all  frequencies,  their  use  is  in  fact  very  often  impractical 
for  high  frequencies,  particularly  for  randomly  excited  complex 
structures.  The  classical  approach  consists  of  determining  the 
natural  modes,  of  calculating  the  responses  of  these  modes  to  a 
specified  excitation  of  interest,  and  of  superposing  these 
responses  to  determine  the  total  structural  response.  Contin¬ 
uous  structures  have  an  infinite  number  of  modes,  but  generally 
only  the  lowest  few  of  these  are  of  importance  in  low  frequency 
vibrations,  so  that  in  these  cases  one  needs  to  consider  only 
those  few  modes.  At  high  frequencies,  however,  a  frequency  band 
of  interest  usually  encompasses  the  resonances  of  a  large  number 
of  modes,  and  one  must  consider  the  responses  of  all  of  these 
modes  in  calculating  the  structural  vibrations  in  that  band. 

Such  classical  multimodal  analyses  may  require  an 
amount  of  computation  which  is  so  great  that  It  exceeds 
the  capacities  of  some  of  today's  largest  computers.  But,  one 
also  encounters  another,  perhaps  more  significant,  limitation 
on  the  utility  of  such  analyses.  Since  the  shapes  of  the 
higher  modes  are  much  more  sensitive  to  details  of  the 
structure  than  are  those  of  the  lower  modes,  one  must  be 
able  to  describe  the  structural  and  material  properties 
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with  greater  precision  in  order  to  be  able  to  perform 
higher  mode  calculations  meaningfully.  The  required 
precision  can  generally  be  achieved  only  for  fictitious 
structures,  in  view  of  the  manufacturing  tolerances  and 
the  variation  in  material  properties  one  inevitably 
encounters  in  realistic  cases. 

If  one  could  perform  a  high-frequency  multimodal 
analysis  meaningfully,  one  would  still  be  faced  with  the 
problem  of  interpreting  the  large  mass  of  results  that  is 
generated  in  such  a  set  of  calculations.  These  results 
would  be  descriptions  of  the  motions,  within  the  frequency 
band  of  interest,  of  a  large  number  (ideally  all)  of  the 
points  of  the  structure.  In  general,  these  motions  will 
be  complex  and  different  for  different  points,  and  one 
will  need  to  perform  additional  calculations,  e.g.  to 
determine  average  values  of  various  response  properties, 
in  order  to  reduce  the  large  amount  of  information  to 
tractable  size. 

The  statistical  energy  analysis  approach  to  structural 
vibrations  problems,  pioneered  by  R.  H.  Lyon,  P.  W.  Smith,  Jr., 
and  I.  Dyer,  was  developed  in  response  to  the  need  for  a 
simple  means  for  understanding  and  estimating  significant 
properties  of  multimodal  vibrations  of  complex  systems. 

This  approach  was  spawned  by  the  realization  that  averaging 
initially,  and  then  carrying  out  calculations  in  terms  of 
average  quantities,  should  lead  to  results  much  more  readily 
than  the  classical  approach,  which  involves  much  initial 
detailed  calculation  and  subsequent  averaging. 

The  statistical  energy  approach  to  some  extent  is 
analogous  to  the  "room  acoustics"  approach  for  dealing  with 
sound  in  architectural  spaces.  In  the  latter  approach  one 
does  not  attempt  to  solve  the  acoustic  wave  equation  in 
detail  for  the  complex  spaces  involved;  rather,  one  analyzes 
the  average  behavior  of  energy  variables  and  interprets  the 
results  in  terms  of  averages  of  dynamic  variables  (e.g., 
pressures,  velocities).  In  the  statistical  energy  approach 
one  similarly  does  not  concern  oneself  with  a  detailed 
solution  of  the  system  equations  of  motion,  but  one  analyzes 
average  energy  distributions  and  relates  those  to  average 
response  variables. 

Although  the  term  "statistical  energy  method"  has  come 
into  somewhat  popular  use  at  this  writing,  there  really 
exists  no  "method"  consisting  of  a  well-defined  technique 
leading  to  a  guaranteed  result  (like  the  Rayleigh-Ritz  or 
the  Holzer  methods,  for  example).  There  does  exist  an 
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approach,  which  -  like  room  acoustics  -  is  one  of  a  set  of 
tools,  the  intelligent  use  of  which  can  give  one  approximate, 
but  simple  and  useful  answers  to  problems  which  could  not 
otherwise  be  solved  realistically. 

The  fundamentals  of  the  statistical  energy  analysis 
approach  are  scattered  in  the  technical  literature,  generally 
are  couched  in  language  which  may  discourage  the  non-specialist, 
and  in  some  cases  have  been  only  implied,  and  not  spelled  out 
in  detail.  Consequently,  this  potentially  highly  useful 
approach  has  not  come  to  be  used  as  widely  as  one  may  have 
expected  on  the  basis  of  its  simplicity  of  application.  It 
is  the  purpose  of  this  report  to  fill  the  apparent  need  for 
a  presentation  of  the  fundamentals  of  the  statistical  energy 
analysis  approach  in  easily  accessible  form,  for  an  explana¬ 
tion  of  its  area  of  applicability  and  for  an  illustration  of 
how  these  fundamentals  may  be  used. 

The  first  of  the  following  sections  reviews  the  concepts 
and  properties  of  modes  of  structural  vibrations,  and  points 
out  some  useful  relations  between  modal  response  and  total 
average  response  properties.  The  second  section  derives  the 
basic  relations  that  govern  the  exchange  of  energy  between 
two  coupled  modes.  The  third  section  generalizes  this  relation 
to  enable  one  to  determine  the  flow  of  energy  from  one  set  of 
modes,  representing  one  structure  or  fluid  system,  to  another 
set  of  modes,  representing  another  such  system.  The  fourth 
section  relates  some  of  the  parameters  introduced  in  the  modal 
analyses  to  more  commonly  employed,  and  more  easily  estimated 
and  measured  system  parameters.  The  fifth  and  final  section 
illustrates  the  application  of  the  statistical  energy  approach 
to  a  special  class  of  problems. 

The  reader  who  does  not  have  the  inclination  to  wade 
through  the  considerable  amount  of  mathematical  manipulation 
which  appears  in  this  report  should  not  be  discouraged.  He 
should  ignore  all  those  details  and  concentrate  on  the  results 
and  conclusions  Indicated  In  the  "Summary  and  Conclusions" 
subsections.  The  strong  point  of  the  statistical  energy 
approach,  after  all,  is  the  simplicity  with  which  the  central 
results  can  be  stated  and  applied. 
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MODES  OF  CONTINUOUS  SYSTEMS 


If  an  undamped  linear  elastic  system  vibrates  in  the 
absence  of  external  forces  in  such  a  way  that  all  points  in 
the  system  move  sinusoidally  in  time  at  the  same  frequency, 
then  the  system  is  said  to  vibrate  in  a  "natural  mode"  -  or 
"mode",  for  short.  This  type  of  vibration  is  a  very  special 
member  of  the  class  of  all  possible  types  of  free  vibrations, 
and  possesses  some  useful  properties,  several  of  which  are 
reviewed  in  the  present  section. 

The  various  points  in  an  undamped  system  vibrating  in 
a  natural  mode  generally  oscillate  at  different  amplitudes, 
but  are  either  completely  in  phase  or  out  of  phase  with  each 
other.  All  points  reach  their  maximum  excursion  (some  in 
the  negative  and  some  in  the  positive  direction)  at  the  ssune 
instant,  all  points  pass  through  their  equilibrium  positions 
simultameously,  aind  all  points  reach  any  fraction  of  their 
maximum  excursions  simultaneously.  Thus,  a  natural  mode 
oscillation  may  be  characterized  by  a  spatial  distribution 
of  deformations  (e.g.,  by  the  distribution  of  deflections  at 
any  suitable  instant),  multiplied  by  a  sinusoid  in  time  and 
a  constant  which  is  indicative  of  the  amplitudes. 

The  spatial  distribution  of  deformations  is  called  a 
"mode  shape" .  Functions  describing  a  given  mode  shape  may 
differ  from  each  other  by  a  multiplicative  constant,  since 
the  definition  of  the  mode  shape  is  arbitrary  to  that  degree. 
However,  once  a  mode  shape  function  has  been  selected,  the 
constant  by  which  one  must  multiply  this  function  so  that  it 
represents  a  given  modal  vibration  is  fully  determined.  This 
latter  constant  is  called  "modal  amplitude";  a  specified  modal 
amplitude  and  mode  shape  function  together  define  a  unique 
distribution  of  deformations. 

The  frequency  at  which  a  system  in  absence  of  external 
forces  vibrates  with  a  given  mode  shape  is  called  the  "natural 
frequency"  corresponding  to  that  mode.  A  continuous  elastic 
system  can  vibrate  in  an  infinite  number  of  different  modes, 
with  an  infinite  number  of  corresponding  natural  frequencies. 

The  discussion  of  modes  presented  here  is  of  necessity 
limited  in  scope  and  rigor.  A  clearly  written  more  complete 
treatment  may  be  found  in  Reference  1;  a  collection  of  infor¬ 
mation  on  the  modal  properties  of  simple  systems  is  available 
in  Reference  2. 
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Mode  Shapes  and  Natural  Frequencies 


The  linear  partial  differential  equations  that  govern 
the  motions  of  undamped  continuous  elastic  svstems  (such  as 
strings,  membranes,  shafts,  beams  and  plates)  in  absence  of 
external  forces  may  be  written  in  the  form 

|x  u  +  K  u  =  0  (l) 


Here  K  is  a  linear  differential  operator  which  involves  only 
the  spatial  coordinates  and  represents  the  stiffness  charac¬ 
teristics  of  the  system,  and  p.  is  a  function  of  only  the 
coordinates  and  describes  the  mass  distribution  of  the  system. 
The  symbol  u  denotes  the  deflection  of  the  system  from  equi¬ 
librium  as  a  function  of  time  t  and  of  the  spatial  coordinates, 
and  the  dot  atop  a  symbol  indicates  differentiation  with 
respect  to  time. 

[For  example,  for  a  shaft  in  torsion  u=u(x,t)  denotes  the 
torsional  angular  deflection  as  a  function  of  the  position 
along  the  shaft  and  of  time.  The  corresponding  operators 
for  Eq.  (l)  are 

M-  =  pJ,  K  =  -  4(ksG  , 

where  p  denotes  the  material  density,  J  the  polar  moment 
of  inertia  of  the  cross-sectional  area,  Ks  the  torsional 
constant  of  the  crosssection,  G  the  shear  modulus  of  the 
material.  In  the  most  general  case  all  of  these  parameters 
may  be  function  of  the  coordinate  x  measured  along  the 
length  of  the  shaft. 

Similarly,  for  a  flat  plate  in  flexure,  u=u(x,y,t)  denotes 
the  displacement  of  the  plate  normal  to  its  midsurface, 
as  a  function  of  coordinates  x,y  along  that  midsurface 
and  of  time.  Here 

H  =  Ps(x,y) ,  K  =  V2(D  V2) 


where  p  denotes  plate  mass  per  unit  surface  area,  and  D 
represents  the  plate  flexural  rigidity,  usually  given  as 

Ehvl2(l-v  )  for  homogeneous  plates  of  thickness  h,  2 
Young's  modulus  E  and  Poisson's  ratio  v.  The  symbol  v 
represents  the  "Laplacian"  operator;  in  the  usual  Cartesian 


coordinates  V  = 


dx‘ 


V 


In  the  most  general  isotropic 


case,  p  and  D  may  both  vary  with  x  and  y. ] 
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Natural  modal  vibrations  are  defined  as  free  (unforced) 
vibrations  in  which  all  points  move  in  unison,  sinusoidally 
in  time,  at  the  same  frequency.  Such  vibrations  correspond 
to  solutions  of  Eq.  (l)  which  are  of  the  form* 


u(x,y,t)  =  ^(x,y)  [B  cos(cot  +  C)]  ,  (2) 

in  which  the  spatial  and  time  dependences  are  separated  and 
where  B  and  C  are  constants.  Substitution  of  Eq.  (2)  into  (l) 
yields  the  time-independent  partial  differential  equation 

~  2 

-  co  [nf/  =  0  (o) 


which  the  function  ^(x,y)  of  the  spatial  coordinates  must 
satisfy. 

One  finds  that  solutions  ifj  of  Eq.  (3),  which  satisfy  the 
boundary  (edge)  conditions  applicable  in  a  given  case  and  which 
are  not  identically  zero,  exist  only  for  certain  values  of  the 
frequency  co.  These  values  are  called  the  natural  frequencies, 
and  the  corresponding  solution  functions  ^(x,y)  are  called  the 
mode  shapes  of  the  system.  To  each  n-dimensional  system 
(n=l,2,3)  there  correspond  n-fold  infinite  sets  of  natural 
frequencies  and  mode  shapes.  To  each  two-dimensional  system 
(plate,  membrane)  there  correspond  doubly  infinite  sets  of 
natural  frequencies  co..  and  mode  shapes  ib (x,y):  1,k  =  1.2.*«* 
which  satisfy  Eq.  (3)Jor  Jk  ’ 

K^Jk(x,y)  -  co^k  mfrJk(x,y)  =  0  (4) 


Orthogonality  and  Normalization  of  Mode  Shapes 

The  mode  shapes  of  systems  whose  edges  are  elastically 
constrained  against  lateral  and  rotational  deformation2  possess 


iFrom  here  on  the  notation  used  is  that  for  two-dimensional 
systems.  Expressions  applicable  to  one-  and  three-dimensional 
systems  can  be  deduced  by  direct  analogy.  (Also,  see  Ref.  2.) 

sZero  and  infinite  values  of  constraint  stiffness  are  permitted, 
but  in  the  general  case  the  boundaries  exert  restoring  forces 
which  are  proportional  to  the  deformations  and/or  restoring 
moments  which  are  proportional  to  the  rotations.  Infinite 
lateral  stiffness  and  zero  rotational  boundary  stiffness 
corresponds  to  the  usual  "simple  support"  conditions.  The 
case  where  both  stiffnesses  are  infinite  corresponds  to  clamped 
edges;  that  xfhere  both  are  zero  corresponds  to  free  edges 
(p.  209,  Ref.  1.)  * 
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the  useful  property  of  "orthogonality".  Mathematically, 
this  means  that  the  mode  shapes  satisfy 


M>(x,y)  tfjk(x,y) 


A 


s 


^j,kl(x,y)  dAs 


nfFJk  for  j=j',  k=k« 

(5) 


0  otherwise 


where  the  integration  is  carried  out  over  the  entire  elastic 
system  (here  shown  as  over  the  entire  area  Ag  of  a  two- 
dimensional  system).  In  other  words,  if  the  integral  of 
Eq.  (5)  involves  two  different  modes,  then  the  result  is 
zero,  hut  if  it  involves  the  same  mode  twice,  then  the 
result  is  a  finite  number.  This  number  has  been  denoted  by 
(m¥jk)>  where  m  denotes  the  total  mass  of  the  system,  and 
satisfies 


^jk  S  II  ^  ^Jk(x,y)  dy 


(6) 


in  accordance  with  Eq.  (5)* 

One  may  observe  that  multiplication  of  ij/  by  an  arbitrary 
constant  does  not  affect  Eq.  (3)  a-nd  that  hence  this  equation 
defines  ^(x.y)  only  within  an  arbitrary  multiplicative  constant, 
Equation  (6)  therefore  involves  a  similar  arbitrariness.  It 
is  generally  useful  for  analysis  purposes  to  normalize  the 
mode  shapes  *1k(x,y),  -  that  is,  to  set  down  some  more  or  less 
arbitrary  rule  for  determining  the  multiplicative  constant  to 
be  used  in  a  given  analysis.  One  commonly  used  rule  (but  by 
no  means  the  only  one)  requires  ¥ jj^l  for  all  j  and  k,  so  that 
nflf..  In  Eq.  (6)  reduces  to  the  total  system  mass  m.  If  this 
rule  is  used,  then  the  mode  shape  functions  satisfy 


If  H(x,y)  *kJ(x,y)  dxdy  =  m  =  ff  \i(x,y)  dx  dy 


(7) 


in  addition  to  Eq.  (4). 
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[For  example ,  for  torsional  vibrations  of  a  uniform  shaft 
clamped  at  one  end  (x=0)  and  free  at  the  other  (x=L) ,  one 
finds 

^.(x)  =V~2  sin[  (2J-1)ttx/2L]  ,  cd^  =  [  (2j-l)ir/2L]VKsG/pJ. 

For  a  uniform  rectangular  plate  extending  from  x=0  to  x=a 
and  from  y=0  to  y=b,  and  simply  supported  on  all  four 
edges  one  finds 

*Jk(x,y)  =  2  sin(  jirx/a)  sin(kiry/b)  , 

=  tt2(  j2/a2  +  k2/b2)  VD/psh  .  ] 

Representation  of  Deformations  in  Terms  of  Modes 

Anv  mathematically  "well-behaved"  system  deflection 
u(x,y,t;  may  be  expressed  as  a  sum 


00  00 

u(x,y,t)  =  Y,  Yj  Ujk^  ^jk^x'y)  (8) 

J=1  k=l 


involving  the  mode  shapes  ,,  and  "modal  amplitude  coefficients" 
U1k(t),  which  depend  only  onKtime  (Refs.  1,3)*  To  evaluate  these 
coefficients  one  may  proceed  in  a  manner  which  is  analogous  to 
that  used  to  obtain  a  Fourier  series  expansion  of  a  given  func¬ 
tion.  If  one  multiplies  both  sides  of  Eq.  (8)  by  M-(x,y)^j  i ki  (x,y) 
and  integrates  over  the  entire  elastic  system,  then  one  finds  in 
view  of  the  orthogonality  relation  (5)  that  only  one  of  the 
integrals  of  the  sum  will  not  vanish.  (Namely,  that  for  which 
j=J*  and  k=k'.)  With  the  normalization  of  Eq.  (7)  the  value  of 
this  non-zero  integral  is  m,  and  thus  one  obtains  the  result 


UJk(t)  =  u(x,y, t )  ^jk(x,y)  M-(x,y)  dx  dy  (9) 


which  permits  one  to  determine  the  coefficient  U.k(t)  which 
corresponds  to  the  mode  shape  ^,,,(x,y)  for  a  given  deformation 
u(x,y,t). 
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Modal  Analysis  of  System  Motions 


The  deflection  of  a  viscously  damped  elastic  system  due 
to  the  action  of  a  distributed  load  p(x,y,t)  is  governed  by  a 
differential  equation  of  the  form 


M-  u  +  y  u  +  K  u  =  p(x,y,t) 


(10) 


where  the  various  terms  have  the  same  meaning  as  in  Eq.  (l), 
and  y  is  an  operator  or  function  which  involves  only  the 
coordinates  and  which  accounts  for  energy  dissipation  (viscous 
damping)  in  the  system. 

In  order  to  determine  the  equations  that  govern  the 
responses  of  the  modes  of  the  system  one  may  substitute  for 
u(x,y,t)  in  Eq.  (lo)  the  series  given  ifl  Eq.  (8).  If  one  uses 
Eq.  (4)  to  replace  the  terms  involving  K  by  others  involving  \i, 
one  obtains 


00 


00 


j=l  k=l 

Now  one  may  again  try  to  use  the  strategem  that  was  employed 
in  the  derivation  of  Eq.  (9).  Namely,  one  may  multiply  both 
sides  of  Eq.  (ll)  by  iw (x,y) ,  integrate  over  the  entire 
system,  and  apply  the  orthogonality  condition.  However,  in 
general  the  damping  operator  y  and  the  mass  distribution  p. 
may  not  be  related  in  any  manner  which  permits  one  to  apply 
to  the  y  terms  the  orthogonality  condition  of  Eq.  (5)  which 
involves  |i. 

Fortunately,  the  special  case  where  y=Pix,  with  p  a  constant, 
turns  out  to  be  a  useful  one,  and  may  serve  as  a  good  approxima¬ 
tion  for  many  kinds  of  lightly  damped  systems.  In  this  case  the 
previously  discussed  strategem  does  work.  Its  use  in  conjunction 
with  the  normalization  of  Eq.  (7)  yields  a  result  which  one  may 
write  as 


(12) 


if  one  defines 
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(13) 


For  a  given  elastic  system  subject  to  a  given  loading, 
the  "modal  mass"  m,  the  "modal  damping  coefficient"  c,  and 
the  "modal  stiffness"  Kjk  are  constant,  whereas  the  "modal 
force"  Fjk  is  a  function  of  time  only.  Equation  (12),  which 
governs  the  modal  coefficients  Ujk,  may  be  recognized  to  be 
the  same  as  the  well-known  equation  that  governs  the  displace¬ 
ment  of  a  linear  single-degree-of-freedom  system  consisting  of 
a  mass  m,  mounted  on  parallel  arrangement  of  a  spring  (with 
stiffness  Kjk)  and  dashpot  (with  viscous  damping  coefficient  c), 
and  subject  to  a  force  Fjk(t)  which  varies  with  time. 

It  is  of  interest  to  note  that  Eq.  (12)  applies  for  each 
mode  (i.e.  for  all  values  of  j  and  k) ,  but  that  the  equation 
for  a  given  j  and  k  does  not  involve  terms  in  other  values  of 
j  and  k.  There  is  then  no  interaction  between  the  modal  motions, 
and  the  modes  are  said  to  be  "uncoupled".  Such  uncoupling  does 
not  occur  in  the  most  general  case,  in  which  orthogonality  cannot 
be  applied  to  the  damping  terms  involving  y,  and  in  which  a  sum¬ 
mation  involving  all  modes  would  appear  in  Eq.  (12)  instead  of 
the  single  c  U.jk  term.  However,  as  mentioned  previously,  one 
may  usually  ignore  this  damping  coupling  for  lightly  damped 
systems . 

Spatial  Averages  of  Deformations;  System  Kinetic  Energy 

If  one  squares  both  sides  of  Eq.  (8),  multiplies  the  result 
by  |i,  and  integrates  over  the  entire  system,  one  obtains 


A 


s 


j=l  k=l  j'=l  k'=l  A 

S 


00 


00 


(14) 


j=l  k=l 
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where  the  second  form  of  the  right-hand  side  follows  from 
orthogonality,  Eq.  (5),  and  from  Eq.  (7).  The  weighted 
spatial  average  value  of  u2,  with  the  averaging  done  over 
the  entire  system  and  with  n  being  used  as  a  weighting 
function,  is  defined  as 


u2(t) 


x,y,t)  n(x,y)  dx  dy  ♦ 


dx 


(15) 


By  comparing  this  expression  with  Eq.  (l4)  and  (7)  one  finds 
that 


00  00 

Zujk(t) 


j=l  k=l 


(16) 


and  thus  that  the  weighted  mean  square  displacement  is  equal 
to  the  sum  of  the  squares  of  all  the  modal  displacements 
Ui^t).  For  uniform  systems,  in  which  the  mass  distribution  u 
is  constant,  the  weighted  mean-square  displacement  is  equal  to 
the  ordinary  mean-square  displacement,  as  may  easily  be  verified 
from  Eq.  (15);  Eq.  (16)  then  applies  to  the  ordinary  mean-square 
displacement. 


If  one  writes  an  expression  like  Eq.  (14)  in  terms  of  the 
velocity  v(x,y,t)  and  the  corresponding  modal  velocity  coeffi¬ 
cients  V^,  then  the  integral  on  the  left-hand  side  of  that 
expression  may  be  seen  to  represent  twice  the  total  kinetic 
energy  T,p  of  the  system.  Thus,  one  may  write 


00  00 


Tip(t )  =  -JJ  v2(x,y,t)  (i  dx  dy  =  ^  ^T^t)  t1?) 


j=l  k=l 


where 


(18) 
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denotes  the  kinetic  energy  associated  with  mode  (j,k). 

Equation  ( 17)  shows  that  the  total  kinetic  energy  is  equal 

to  the  sum  of  all  of  the  modal  kinetic  energies ,  and  also  to 

12  2 

m  v  ,  where  the  weighted  mean-square  velocity  v  is  defined 

analogously  to  u  in  Eq.  (15). 

Expressions  for  the  time-  and  space-  averages  may  be 
obtained,  of  course,  simply  by  averaging  both  sides  of 
Eqs.  (l6) -(l8)  with  respect  to  time. 


Summary  and  Conclusions 

The  definition  of  modes  of  elastic  systems  and  the 
properties  of  such  modes  have  been  reviewed  in  rather  general 
terms  in  the  present  section. 

It  has  been  indicated  that  any  physically  reasonable 
displacement,  velocity,  etc.,  distribution  can  be  represented 
as  a  superposition  of  corresponding  modal  distributions. 

Orthogonality  of  the  modes,  which  applies  for  systems  with 
elastically  constrained  boundaries  (and  thus  for  all  the  usually 
considered  ideal  cases  and  for  most  of  those  generally  of  prac¬ 
tical  interest),  has  been  shown  to  lead  to  uncoupled  equations 
of  motions  for  the  individual  modes,  except  where  interaction 
occurs  due  to  damping. 

The  uncoupled  modal  equations  of  motion  have  been 
demonstrated  to  be  precisely  of  the  same  form  as  those  for 
independent  single-degree-of -freedom  mass-spring-dashpot 
systems.  Thus,  by  use  of  the  concept  of  modes  one  may  replace 
the  complex  problem  of  determining  the  motions  of  a  continuous 
distributed  elastic  system  by  the  much  simpler  problem  of  eval¬ 
uating  the  motions  of  single-degree-of-freedom  systems  which 
represent  the  modes  of  the  elastic  system.  The  motions  of  the 
total  system  may  then  be  obtained  by  superposition  of  all  of  the 
modal  motions,  spatial  averages  of  motion  quantities  may  easily 
be  determined  from  the  corresponding  modal  amplitude  coefficients, 
and  the  total  kinetic  energy  of  the  system  may  be  found  simply  by 
adding  all  of  the  modal  kinetic  energies. 
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POWER  FLOW  BETWEEN  TWO  MODES 

In  many  practical  problems  involving  the  vibrations  of 
complex  elastic  systems,  particularly  where  the  excitations 
and  responses  are  random,  one  is  more  concerned  with  under¬ 
standing  the  interactions  of  various  parts  of  complex  elastic 
systems  and  with  estimating  average  vibration  levels,  than 
with  obtaining  more  detailed  information  about  the  system 
vibrations.  One  finds  that  energy  analyses  are  particularly 
useful  in  such  cases.  These  analyses  tend  to  be  simple,  since 
they  can  make  use  of  the  easily  applied  concept  of  conservation 
of  energy.  And,  as  evident  from  the  previous  section,  if  one 
knows  the  kinetic  energy  of  a  system  one  knows  its  mean-square 
velocity. 

Before  concerning  oneself  with  the  interaction  of  two 
(or  more)  interconnected  elastic  systems,  where  each  system 
contains  many  modes,  one  may  do  well  to  develop  an  understanding 
of  the  interaction  of  Just  two  modes.  It  is  clear  that  in  the 
steady  state  the  vibrational  energy  of  a  mode  of  a  two-mode 
system  is  affected  by  the  power  supplied  to  the  mode  from 
external  sources,  by  the  power  dissipated  by  the  mode,  and  also 
by  the  power  that  the  mode  in  question  receives  from  (or  supplies 
to)  the  other  mode.  The  directly  supplied  and  the  dissipated 
powers  can  generally  be  evaluated  and  related  to  modal  parameters 
and  energies  with  relative  ease.  In  order  to  make  an  energy 
approach  useful,  however,  one  also  needs  expressions  which  relate 
the  mode-to-mode  power  flow  to  modal  energies.  It  is  the  purpose 
of  the  present  section  to  develop  such  expressions. 

It  has  been  pointed  out  that  a  mode  of  an  elastic  system 
behaves  precisely  like  an  ideal  single-degree-of-freedom  system 
consisting  of  the  modal  mass  mounted  on  a  spring  (having  the 
modal  stiffness)  and  a  dashpot  (with  the  modal  damping  coeffi¬ 
cient).  Consequently,  the  energy  transfer  between  two  modes  is 
studied  here  by  analyzing  the  corresponding  behavior  of  two 
coupled  mass-spring-dashpot  systems,  as  shown  schematically  in 
Fig.  1.  In  order  to  facilitate  the  analysis,  and  since  linear 
systems  are  of  primary  interest,  all  of  the  present  discussion 
is  restricted  to  linear  coupling. 

The  concepts  of  impedance  and  admittance  are  reviewed  in 
the  first  of  the  following  sections.  Thereafter  the  restrictions 
which  the  assumption  of  conservative  coupling  places  on  these 
quantities  are  delineated.  (Conservative  coupling  is  coupling 
which  neither  dissipates  nor  supplies  energy,  and  thus  is  of 
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primary  interest  where  one  wishes  to  study  the  energy  flow 
between  systems.)  Next  are  developed  expressions  for  the 
modal  kinetic  energies.  Finally,  it  is  demonstrated  that 
the  time-average  power  flow  between  two  modes  is  proportional 
to  the  difference  in  the  time-average  modal  kinetic  energies 
under  some  conditions  which  are  often  approximated  in  practice. 


Admittances  and  Impedances 

In  order  to  maintain  as  much  generality  as  possible,  it 
is  useful  initially  to  place  no  restriction  on  the  coupling 
element  L  of  Fig.  1,  except  for  requiring  it  to  be  linear. 

This  restriction  implies  that  L  can  give  rise  only  to  linear 
terms  in  the  equations  of  motion,  so  that  response  variables 
(displacements,  velocities,  etc.)  are  proportional  to  excita¬ 
tion  variables  (forces,  rate  of  change  of  force,  etc.), 
regardless  of  any  of  the  amplitudes  involved. 

The  properties  of  linear  systems  subject  to  purely  sinu¬ 
soidal  excitations  may  be  studied  conveniently  by  use  of  phasor 
notation  (Ref.  4).  In  this  notation  each  sinusoidally  varying 
dynamic  variable  is  represented  by  a  complex  quantity  or  "phasor", 
and  the  actual  time  functions  are  obtained  by  taking  the  real 
part  of  the  product  of  the  corresponding  phasor  and  e-^*,  where 
co  denotes  the  radian  frequency,  t  time,  and  i  =  n/~-1.  Thus,  for 
example,  if  V-t  is  the  phasor  corresponding  to  the  velocity  V.(t), 
then  0  3 

Vj(t)  =  Re[V^  e"lu)t]  =  Re[V^]-cos  cot  +  ImtV^-sin  cot  .  (19) 

In  view  of  the  linearity  of  the  two-mass  system  of  Fig.  1 
the  velocity  phasors  of  the  two  masses  obey  the  relation 


=  Y 

=  Y 


11 

+ 

Y12 

F 

^-2 

21 

*1 

+ 

Y22 

F 

^2 

(20) 


where  F,  and  F^  denote  phasors  corresponding  to  the  forces  F, 
and  F  The  quantities  Yrs  are  called  "admittances";  from 
Eqs.  ^(20)  it  is  evident  that  Yni  =  Vj/Fj  if  Fp  =  0,  Y10  =  V -1  /Fo 
if  Fj_  =  0,  Y21  =  if  Fg  =  0,  etc. 
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By  solving  Eqs.  (20)  for 
may  write 


11  2i 


V 

21  -1 


+  Z 
+  Z 


12  —2 
V 

22  -2 


and  F2  one  finds  that  one 


(21) 


where 


Z11~Y22//D*  Z22  Yll//D,  Z12  ~Y12/^>’ 

^Y11Y22  "  Y12Y21 


Z2r"Y21^ 


(22) 


The  Zrs  quantities  are  called  "impedances".  From  Eqs.  (21) 
one  may  observe  that  Zji=F if  V2=0,  Z].2=.El/^2  Yn  =0,  etc. 
Alternately,  by  solving  Eqs.  (2l)  for  and  V2  one  obtains 
Eqs.  (20)  and  finds  that  the  admittances  may  be  expressed  in 
terms  of  the  impedances  as 


yH=Z22/S,  Y22=Zn/E,  Y12=-Z12/E,  Y21=-Z21/E 


(23) 


E"Z11Z22"Z12Z21_1/1)  * 


The  impedances  and  admittances  in  general  are  complex 
quantities  and  functions  of  the  frequency  od.  They  depend  in 
general  on  the  uncoupled  system  parameters,  as  well  as  on  the 
coupling  element. 

Implications  of  Conservative  Coupling 

In  order  to  study  the  flow  of  energy  between  systems  it 
is  useful  to  introduce  the  restriction  of  conservative  coupling 
in  addition  to  that  of  linear  coupling;  that  is,  to  postulate 
coupling  which  neither  dissipates  nor  supplies  energy.  Only 
for  such  coupling  is  the  energy  leaving  one  of  two  coupled 
systems  equal  to  the  energy  reaching  the  other,  and  only  for 
such  coupling  can  one  speak  unambiguously  of  energy  flow  from 
one  system  to  the  other.  The  present  section  consequently  is 
concerned  with  the  restrictions  which  are  imposed  on  the  various 
admittances  and  impedances  by  the  requirement  that  the  coupling 
be  conservative. 
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Conservation  of  energy  demands  that  in  the  steady  state 
the  total  work  done  per  cycle  by  the  external  forces  be  Just 
equal  to  the  total  energy  dissipated  per  cycle.  In  view  of 
the  well-known  relations  between  power  and  forces  and  veloci¬ 
ties  in  complex  notation,  as  also  summarized  in  Ref.  5,  this 
statement  of  energy  conservation  may  be  written  as 

|  Ret^V*]  +  |  RelPgVJ]  =  flvj2  Re[Zu]  +  i^l2  Re[Z22)  , 

(24) 


where  the  asterisk  indicates  the  complex  conjugate  of  the 
quantity  to  which  it  is  appended. 


Substitution  for  F-,  and  Fp  from  Eqs.  (21),  replacement 
I  Yi  I  ^  by  YlYl>  and  elimination  of  terms  common  to  both 
sides  of  the  resulting  equation,  permits  one  to  reduce  the 
foregoing  relation  to 


Re[Z12  'k  n  +  Z21  2i  S|J  -  O  .  (25) 

Since  Vp— l  =  (Y^Y^)**  this  result  implies 
(Re[Z21]  +  Re[Z12])  ReQ^V*)  - 

(26) 

(lm[Z21]  -  Im[Z12J)  ImQ^V*)  =  0. 


But  this  relation  must  hold  for  all  V-,  and  Vp,  and  thus  for 
all  values  of  Re(V]V*>)  and  Im(VjV£).  Consequently  the  coeffi¬ 
cients  of  these  two  quantities  in  the  foregoing  expression 
must  vanish,  and 


Re[Z21]  =  -  Re[Z12] 
Im[Z21]  =  Im[Z  12] 


(27) 


Conservative  coupling  thus  implies  that  the  "transfer 
impedances"  satisfy 


Z 


12 


Z* 

*21 


(28) 
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and  in  view  of  Eqs.  (22)  or  (23),  that  the  transfer  admittances 
obey 


(29) 


Time-Average  Modal  Kinetic  Energies 

Since  it  is  desired  to  develop  relations  which  express  the 
power  flow  between  modes  in  terms  of  the  time-average  modal 
kinetic  energies,  it  is  of  interest  first  to  study  these  energies 
from  a  rather  general  viewpoint. 

Clearly,  the  time-average  value  of  the  kinetic  energy  Tj  of 
a  mode  j  is  given  by 

<Tj>  =  <V®(t:)>  (30) 

where  m-.  denotes  the  modal  mass  (or  the  mass  of  the  system 
of  which  j  is  a  mode)  and  Vj  the  modal  velocity,  and  where 
the  brackets  <•••>  indicate  averaging  with  respect  to  time. 

The  time-average  value  of  the  square  of  the  modal  velocity 
is  also  known  as  the  mean-square  modal  velocity;  it  is  inde¬ 
pendent  of  time,  in  the  steady  state.  The  problem  of  deter¬ 
mining  the  average  kinetic  energies  thus  corresponds  to  that 
of  finding  the  mean-square  velocities,  and  these  may  be  found 
from  the  velocities  per  se. 

If  the  velocity  of  mass  1  of  Fig.  1  due  to  a  unit  impulse 
force  5 ( t )  acting  on  this  mass  is  given  by  y11(t),  then  the 
velocity  V,(t)  due  to  a  general  force  F1(t)  acting  on  mass  1 
may  be  expressed  as  a  convolution  (Refs.  1,6); 


Vt)  =  /  yil(ct)  F1(t-<x)  .  (31) 


Here  a  is  a  dummy  time  variable,  which  disappears  after  the 
integration  is  performed  and  limits  are  substituted.  Similarly, 
if  the  velocity  of  mass  1  due  to  a  unit  impulse  force  acting  on 
mass  2  is  given  by  y^t),  then  the  velocity  of  mass  1  due  to  a 
general  force  F2(t)  acting  on  mass  2  may  be  written  as 

00 

Vt)  =  f  y12(P)  Pg(t-P)  dP  ,  (32) 


where  P  is  again  a  dummy  variable. 
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According  to  the  superposition  principle,  the  response 
produced  by  a  combination  of  excitations  acting  on  a  linear 
system  is  equal  to  the  sum  of  the  responses  due  to  the  indi¬ 
vidual  excitations.  Thus,  the  velocity  Vj(t)  of  mass  1  due 
to  the  combination  of  a  force  F-^t)  acting  on  mass  1  and  a 
force  Fg(t)  acting  on  mass  2  may  be  expressed  as 

00  00 

Vi(t)  =  f  yxl(a)  Px(t-a)  d a  +f  y12(P)  Pg(t-P)  dP  .  (33) 


The  mean-square  value  of  V-^(t)  is  obtained  by  squaring 
the  foregoing  equation  and  averaging  over  all  time.3  If  the 
two  integrals  of  Eq.  (33)  are  denoted  by  J11=J,1(t)  and 
Jl2=Jl2(t),  respectively,  one  may  write 

<V^>  =  <(J11+J12)2>  =  <J21>  +  <J22>  +  2<J±1  J12>  .  (34) 


The  mixed  product  term  may  be  expressed  as 
00  00 

<J11  J12>  =  II  yll(o0  y12(P)  <F1(t-a)  P2(t-P)>  da  dp  ,  (35) 

—  00  —00 


where  the  averaging  brackets  have  been  placed  only  around  those 
quantities  which  involve  the  variable  t,  with  respect  to  which 
the  averaging  is  to  be  performed.  If  F-.  and  Pp  are  uncorrelated, 
that  is  if 

00 

f  F1(t-a)  F2(t-fS)  <Jt  =  0  (36) 

-00 

for  all  a,  p,  then  <P.  (t-a)  F  (t-P)>  =  0,  and  consequently 
<J11  J12>  -  °- 


3In  most  cases  one  obtains  an  adequate  approximation  by 
averaging  over  many  periods  of  the  lowest  frequency  contained 
in  V-Jt). 
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The  first  term  of  Eq.  (34)  may  be  written  analogously 
to  Eq.  (35)  as 

00  00 

<J^>  =//  yn(a)  yn(P)  <Fx(t-a)  P1(t-P)>  da  dp  .  (37) 


Here  the  term  within  the  averaging  brackets  may  be  recognized 
as  the  autocorrelation  function  R^  of  P^(t);  this  autocorrela¬ 
tion  function  obeys  the  relation  (Ref.  6,  pp.  20,  32-34) 

Rx(a-p)  =  <P1(t-a)  P1(t-P)>  =  <P1(t)  P1(t-a+p)>  = 

00 

4/  S1(03)eia)(a_p)  dco  (38) 


where  is  defined  as  the  mean  square  spectral  density  of  the 
stationary  random  force  F^(t)  and  is  related  to  the  autocorrela¬ 
tion  function  as 


Si(co)  =  f  R1(r)e"iayr  dr  .  (39) 


The  functions  Sn(to)  and  R1(r)  are  a  Fourier  transform  pair. 
So  also  are  the  impulse  response  and  frequency  response  functions 
y^k(t)  and  Yjk(co),  which  for  all  j,k,  obey* 


4It  should  be  pointed  out  that  several  different  valid  Fourier 
transform  definitions  may  be  used.  These  differ  only  in  the 
coefficients  attached  to  the  integrals,  but  all  result  in  a 
factor  of  1/2tt  for  a  complete  transformation  cycle.  E.g.,  some 
authors  use  l/Vl^ir  in  each  transformation.  Here  1/2tt  has  been 
used  where  the  integrations  are  performed  with  respect  to  o 
and  unity  where  the  integrations  are  performed  with  respect  to  r. 

It  is  also  important  to  note  that  Eqs.  (38)-(40)  involve  both 
positive  and  negative  frequencies  (and  times),  whereas  only 
positive  frequencies  have  physical  meaning.  An  experimentally 
determined  spectrum  (involving  only  positive  co)  may  be  converted 
to  a  spectrum  corresponding  to  Eqs.  (38)  and  (39)  by  halving 
the  measured  spectrum  at  each  oo  and  plotting  the  result  at  both 
positive  and  negative  numerical  values  of  cu  (Ref.  12,  pp.  11,12). 
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00 

yjk(t)  =^/V“)elu>t  > 

-00 


00 

YJk(0))  =  /* jk(T)e~i<OT  dT 

-00 


(40) 

00 

YJk(“>  -  /  rJk<T)ei‘w'  dT- 
-00 


Substitution  of  Eq.  (38)  into  (37)  gives 
00  00  00 

<J^>  =/  /  /  yn(a)  Yn^)  si(“)  eico(a_P)  da>  da  dP 

-00  —00  -00 

00  00  00 

=  fsiM[  J  y11(a)  eia:0i  da  J  yn0)  e_lcuP  dpj  dco 

—  00  —  00  —  00 

(41) 

00  00 

=  J  S1(co)  Y^co)  Y11(a>)  doo  =  j  S-^co)  |  Y^iia)  |  2  dco  . 

-00  -00 


Analogously  one  may  write 

00 

<J12>  =  f  I  Yi2(co)  1 2  dco 

-00 


(42) 


where  Sp(co)  denotes  the  mean-square  spectral  density  of  the 
force  Fp( t) . 

Substitution  of  the  foregoing  results  into  Eq.  (34) 
permits  one  to  write  the  mean-square  velocity  <V^>  of  mass  1  as 

00  00 

<V}>  =  f  S1(co)|Y11(co)|2  dco  +  /  S2(co)  |yi2|2  dco  ;  (43) 

-00  —00 
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and  to  set  down  by  analogy  the  mean-square  velocity  of  mass  2  as 


00  00 

<V2>  =/S  i(<^)  I  y2i  I  2  dco  +  /  s2{o>)  |y22[2  dco 


-00 


(44) 


Calculation  of  Net  Power  Flow  Between  Modes 

by  Superposition  of  Plows 

In  order  to  study  the  power  flow  (i.e.,  the  rate  of  flow  of 
energy)  from  one  mode  to  another  it  is  convenient  to  consider 
one  force  to  act  at  a  time,  and  then  to  add  the  power  flows 
obtained  in  the  two  cases  to  obtain  the  net  power  flow.  Such 
addition  of  power  flows  is  permissible  for  statistically  inde¬ 
pendent  forces  F-^  and  F2,  as  demonstrated  below. 

Since  the  system  under  consideration  is  linear,  the  force  F 
exerted  on  mass  1  by  the  coupling  element  (or  by  the  mass  on  the0 
coupling  element)  obeys 

00  00 

Pc(t>  =/h ll<°o>  dao  +  /h12<Po>  yt-V  aPo  (“S) 


in  analogy  with  Eq.  (33).  Here  h-,-,  and  h,p  represent  impulse 
response  functions;  h11(t)  denotesxthe  coupling  force  F (t) 
resulting  from  a  unit  impulse  of  force  acting  on  mass  1,  h12(t) 
similarly  denotes  the  F (t)  resulting  from  such  an  impulse 
acting  on  mass  2,  a  ana  0  are  dummy  time  variables. 


The  time-average  power  P12  supplied  to  the  coupling  element 
may  be  obtained  by  multiplying  the  instantaneous  force  as  given 
by  Eq.  (45)  by  the  instantaneous  velocity  as  given  by  Eq.  (33), 
and  averaging  with  respect  to  time.  One  obtains  a  result  which 
one  may  write 


P12  "  <'o<‘>  Vl<t»  "  P12a  -  P21b  + 


AP 


(46) 


where 


00  00 

P12a  -/  /hll<“ o>  yll(a>  <Pl(t-“o)  Pi(t-°)>  d“o  da 
—  00  -00 


21b 


h12(Po)  y12(P)  <F2(t-Po)  P2<t-e)>  dPo  dP 


—  00  -00 


(47) 


00  00 

AP  -/  /  hn<00>  yi2<p)  cMt-V  p2<t-‘3»  dfJ 


00  00 

■/  /h12(P0>  yll<“>  <'2(^0)  Pl(t-a»  <#><>  d“ 


For  conservative  coupling,  the  time-average  power  flowing 
into  the  coupling  element  from  mode  1  at  steady  state  is  exactly 
equal  to  the  power  flowing  into  mode  2,  so  that  f-^2  of  Eq.  (46) 
represents  also  the  power  flow  from  mode  1  to  mode  2.  If  F1(t) 
and  ^(t)  are  uncorrelated,  so  as  to  satisfy  Eq.  (36),  then  the 
time  averages  appearing  in  the  integrands  of  the  foregoing 
expression  for  AP  vanish  and  AP=0.  Since  ?12a  depends  only  on 
the  excitation  F^(t),  it  may  be  interpreted  as  the  power  flow 
from  mode  1  to  mode  2  due  to  the  action  of  Fi(t)  by  itself 
(i.e.,  with  F2=0).  P2lb  be  similarly  interpreted,  except 
that  the  minus  sign  included  In  its  definition  changes  the 
direction  considered  that  of  positive  power  flow;  P2ib  thus 
represents  the  power  flow  from  mode  2  to  mode  1  due  to  the 
action  of  F2(t)  by  itself. 

In  view  of  the  foregoing  discussion  one  may  thus  write  the 
net  time-average  (steady-state)  power  flow  from  mode  1  to  mode  2 
as 


P  =  P  -  P 
r12  12a  21b 


(48) 


where  now  Pipa  denotes  the  flow  from  mode  1  to  mode  2  if  only 
mode  1  is  excited,  P2Th  denotes  the  flow  from  mode  2  to  mode  1 
if  only  mode  2  is  excited,  and  where  the  two  excitations  are 
taken  to  be  uncorrelated. 
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One  may  calculate  the  power  flow  between  modes  directly 
by  carrying  out  the  integrations  indicated  in  Eqs.  (47),  as 
Scharton  has  in  essence  done  for  a  particular  kind  of  coupling 
(Ref.  7).  However,  one  may  retain  a  greater  degree  of  gener¬ 
ality  and  present  a  discussion  which  is  relatively  easily 
understood  by  taking  an  approach  based  on  Eq.  (48j  and  on 
energy  balance  arguments,  as  in  the  following  paragraphs. 

General  Expression  for  Power  Flow 

If  only  a  force  F^  is  present  (that  is,  if  F2=0),  then 
mass  2  is  excited  only  through  the  coupling  element,  and  in 
the  steady  state  all  of  the  power  dissipated  by  the  damping 
of  system  2  must  be  supplied  to  this  system  from  system  1  via 
the  coupling  element.  An  analogous  statement  applies  if  F1=0 
and  F2^0.  Let  the  time-average  power  flow  from  system  J  to 
system  k  be  denoted  by  and  let  the  added  subscript  "a" 

designate  quantities  corresponding  to  the  F2=0  case,  and 
subscript  ,rb"  to  the  Fn=0  case.  In  view  of  the  previous 
energy  balance  arguments  and  of  Eqs.  (43)  and  (44)  one  may 
write5 


P12a  "  C2<V2a>  “  c2 


00 

/  S1(<o)  |Y21|2  dco 


-00 

00 


P21b  -  cl<V?b>  -  C1  /S2<“>  1*12 1 2 


do 


(49) 


where  S,  and  S?  denote  the  mean  square  spectral  densities  of 
the  forces  F-,  and  F^.  In  view  of  Eqs.  (48)  and  (49),  the  net 
power  flow  from  system  1  to  2,  when  both  external  forces  are 
acting,  obeys  the  equation 


(50) 


5Note  that  Cj  =  Re[Zi:L]for  the  system  of  Fig.  1,  if  the  coupling 
Is  conservative.  In  view  of  the  definition  of  Z,,,  as  indicated 
following  Eq.  (22),  Z is  measured  with  V2=0;  hence  no  (dissi¬ 
pative)  elements  in  system  2  can  contribute  terms  to  Z-q.  Since 
the  coupling  is  conservative  it  cannot  contribute  real  terms  to 
Zjj.  Analogously,  c2  =  Re[Z22]. 
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Relation  Between  Power  Flow  and  Kinetic  Energies 


If  Fn  and  F2  are  "white  noise"  forces,  that  is  if 
S^cd)  =  st  and  S2 (<^>)  =  S2  are  constant  (independent  of 
frequency),  then  one  may  introduce  the  notation 


(51) 


and  rewrite  Eqs.  (33)  and  (3^)  as 


<7^  -  S1  In  +  S2  I12 


2>  —  S1  I2j  +  Sg  ^22 


and  Eq.  (50)  as 


(52) 


P12  =  I12^C2S1"C1S2^  *  (53) 

By  solving  Eqs.  (32)  for  and  Sg  in  terms  of  <VP>  and  <v|> 

and  introducing  the  result  into  Eq.  (53) >  one  obtains  an  expres¬ 
sion  for  the  power  flow  P12  in  terms  of  the  mean  square  velocities. 
By  expressing  the  mean-square  velocities  <Y?>  in  terms  of  the 
time-average  kinetic  energies  <T.>  as  indicated  in  Eq.  (30), 
one  finds  that  one  may  then  writeJ 


P12  =  '  021<T2>  (54) 

where 

21  &l  +  I  ) 

12\m,  x22  m,  x12/ 

<*>12  E  - 1 - T1 - 

Xcl  T  T 

■*■11  22  X12 

(55) 

2  hJk  Jn +  i  Jiz) 

d)  =  _ £ _ £ _ 

^21  2 
T  T  T  C~ 

^11  x22  ~  x2 
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and  where  one  has  made  use  of  the  fact  that  I,p=I21,  as 
implied  by  Eq.  (29). 

Power  Plow  for  Weak  Coupling 

The  two  systems  of  Fig.  1  may  be  said  to  be  weakly 
coupled  if  the  transfer  admittance  integrals  I22=I21  have 
values  which  are  much  smaller  than  those  of  the  "self- 
admittance"  integrals  and  122'  and  the  transfer 
impedances  and  admittances  z21,  Y12>  Y21>  do  not 

affect  the  values  of  122  and  I22  significantly. 


Then,  substitution  of  the  first  of  Eqs.  (23)  into 
Eq.  (51)  yields 


'22 


Z11Z22-Z12Z21 


dOi 

(56) 


+  C1  +  k^ioo)-1 


clml 


Note  that  the  approximations  indicated  in  Eq.  (56)  apply 
even  if  the  transfer  impedances  and  admittances  are  not 
strictly  negligible  (as  compared  to  the  self  terms)  over 
the  entire  frequency  range,  provided  that  the  value  of  the 
integral  is  not  significantly  affected  by  contributions 
corresponding  to  transfer  term  values  which  are  not  negli¬ 
gible.  For  example,  if  Z22  and/or  Z22  are  small  compared 
to  Z22  and/or  Z22  in  certain  frequency  ranges  (due  to 
resonances  of  masses  1  and/or  2) ,  then  the  approximations 
of  Eq.  (56)  still  hold,  provided  that  these  frequency 
ranges  are  narrow  enough. 

By  use  of  the  previously  stated  requirement 
I12  «  In»  122  for  loose  coupling  and  of  Eq.  (56)  one  finds 
that  one  may  restate  Eqs.  (54)  and  (55)  as 


^12  “  ^21  ~  2 


tt  C1C2I12  ^c 


(57) 


P12  =  V<T1>  -  <T2>)  * 
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Power  Flow  for  Stiffness  and  Gyroscopic  Coupling 

If  the  coupling  element  consists  of  an  ideal  stiffness 
(represented  by  a  spring  constant  k  )  and  an  ideal  gyroscopic 
element6  (represented  by  a  force/velocity  coefficient  c„) ,  one 
finds  that 


Z 

Z 

Z 

Z 


11 

22 

12 

21 


m. 


m. 


g 


(ico)  +  ^  +  (k1+kc)(i(u) 
(io>)  +  c2  +  (k2+kc)(ioj) 
-  kgdco)-1 

cg  -  V1")'1  • 


-1 

-1 


(58) 


One  may  substitute  these  impedances  into  Eqs.  (23)  to  find  the 
various  corresponding  admittances,  and  one  may  then  evaluate 
the  integrals  In,  1 22 >  ^12  ex&ctly  (by  use  of  the  theorem  of 
residues,  as  summarized  in  Ref.  9>  for  example,  or  from  the 
tabulated  expression  given  on  p.  72  of  Ref.  6).  After  substi¬ 
tution  of  the  results  into  Eqs.  (55)  and  a  considerable  amount 
of  algebraic  manipulation  one  finds  that  the  two  coefficients 
0,2  and  021  are  here  exactly  equal,  regardless  of  the  magnitude 
of  the  coupling,  and  are  given  by 


^12  "  *21 


-♦e  - 


m,m 


?  2 

k  a  +  c  a 
c  k  g  c 

2\2 


lm2  Vc  +  «-<) 


2  1' 


(59) 


6A  gyroscopic  element  is  defined  as  one  which  produces  a 
negative  force  on  mass  2  due  to  a  positive  velocity  of 
mass  1,  if  it  results  in  a  positive  force  on  mass  1  due 
to  a  positive  velocity  of  mass  2.  (See  Ref.  8.)  Note 
that  a  viscous  damping  element  results  in  forces  of  the 
same  sign  on  both  masses  due  to  positive  velocities,  and 
that  gyroscopic  coupling  is  conservative  since  the  Z,p  and 
Z21  Eqs*  (58)  satisfy  the  requirement  of  Eq.  (28). 
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where 


a 


k 


a 

c 


2(^0^ 


c^k2  “f*  c^k-^ 

(ki+kc)/mi  , 
2  m-jO^  , 


"t  £  2^2  ^ 

a)lCJ°2  (  ^  1°°2  +  ^  2tul  ^ 
2 

o2  =  (k2+kc)/m2 
^2  =  2  m2co2 


(60) 


The  definitions  of  Eqs.  (6o)  were  chosen  so  that  co.,  represents 
the  natural  frequency  of  the  j^h  mass  (or  mode)  if J the  other 
mass  is  kept  from  moving,  and  so  that  £  «  denotes  the  ratio  of 
the  damping  coefficient  c*  to  the  critical  damping  coefficient 
2  nijCDj  of  the  jt*1  mass.  J 

The  foregoing  result  is  found  to  agree  with  that  which 
Scharton  (Ref.  7)  derived  for  pure  spring  coupling  by  calcu¬ 
lating  explicitly  the  power  flow  through  the  coupling  element. 
For  small  coupling  the  foregoing  result  also  reduces  to  that 
which  was  obtained  by  Lyon  and  Maidanik  (Ref.  10)  by  means  of 
an  explicit  power  flow  calculation  based  on  small  coupling 
approximations . 

Finite  Frequency  Band  Excitation 

The  previous  discussion  was  based  on  force  spectra  S-. 
and  So  that  are  assumed  constant  for  all  frequencies.  The 
results  obtained  there  similarly  apply  to  power  flow  and 
system  energies  involving  the  entire  infinite  frequency  band. 
However,  a  little  reflection  permits  one  to  interpret  these 
results  also  for  cases  where  only  a  finite  frequency  band  is 
of  interest. 


One  may  obtain  expressions  pertaining  to  only  a  finite 
frequency  band  by  replacing  the  infinite  limits  appearing  in 
the  integrals  of  Eqs.  (4Q)-(5l)  by  the  frequencies  limiting 
the  band.  Equations  (52)-(55)  then  still  apply,  although  the 
various  integrals  that  enter  them  now  are  different. 
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The  magnitude  of  the  admittance  Y(o>)  of  a  not  too  highly- 
damped  single-degree-of-freedom  system  is  much  greater  near 
the  system  natural  frequency  (resonance)  than  at  other  fre¬ 
quencies,  so  thatoothe  greatest  contribution  to  the  value  of 
the  integral  I  =_/|Y(a>)|1 2  dcu  is  due  to  the  frequency  region 
near  the  system  resonance.  Consequently,  the  value  of  the 
integral  will  be  changed  only  slightly  if  the  infinite  limits 
are  replaced  by  finite  ones,  provided  that  the  interval  of 
integration  (i.e.,  the  frequency  band  under  consideration) 
still  encompasses  the  system  resonance.  (See,  for  example. 

Ref.  11,  or  pp.  21ff ,  29ff  of  Ref.  5.) 

Similar  reasoning  applied  to  a  two-degree-of-freedom 
system  (made  up  of  two  coupled  single-degree-of-freedom 
systems),  which  has  two  resonances,  leads  one  to  conclude 
that  the  results  of  Eqs.  (56),  (57),  (59)  and  (60)  also  are 
good  approximations  if  the  excitation  frequency  bands  are 
finite,  but  encompass  both  system  resonances. 

In  view  of  the  previously  discussed  predominance  of  the 
responses  near  the  resonances,  the  aforementioned  results 
may  also  be  taken  to  be  good  approximations  for  cases  where 
the  excitation  spectra  are  not  flat.  All  that  is  required 
for  these  approximations  to  hold  is  that  the  excitation 
spectra  exhibit  no  violent  peaks;  i.e..  no  peaks  that  are  of 
the  same  order  as  the  peaks  in  the  |y|2  curves. 

Summary  and  Conclusions 

It  has  been  shown  in  this  section  that  the  steady-state 
time-average  power  flow  from  one  mode  to  another  is  propor¬ 
tional  to  the  difference  between  the  time-average  kinetic 
energies  of  the  two  modes,  provided  that  both 

1)  the  coupling  between  the  two  modes  is 

-  linear  (giving  rise  to  a  linear  differential  equation), 

-  conservative  (neither  supplying  nor  extracting 
mechanical  energy),  and 

-  light  and/or  purely  spring-like  and/or  gyroscopic; 

and 

2)  the  forces  acting  on  the  two  modes  are  uncorrelated 
and  have  spectra  that  are  relatively  flat  (as  com¬ 
pared  to  the  system  admittance  spectra)  within  the 
frequency  band  encompassed  by  the  resonances  of  the 
coupled  system. 
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Expressions  which  give  the  constant  of  proportionality 
<t>  between  power  flow  and  kinetic  energy  difference  in  terms 
of  the  characteristics  of  the  modes  and  of  the  coupling 
elements  have  been  derived.  These  expressions  are  given  in 
Eqs .  (57)  and  (59). 

Since  the  square  of  the  difference  between  the  squares 
of  the  modal  resonance  frequencies  appears  in  the  denominator 
of  Eq.  (59),  one  may  note  that  the  proportionality  constant  0C 
decreases  rapidly  as  the  frequency  difference  increases.  Thus, 
the  power  flow  between  two  modes  whose  resonance  frequencies 
are  nearly  the  same  will  be  much  greater  than  the  flow  between 
similar  modes  whose  resonances  are  widely  separated. 

The  relation  between  power  flow  and  kinetic  energies, 
which  has  been  established  in  this  section,  permits  one  to 
analyze  the  average  dynamic  behavior  of  coupled  modes  rather 
simply  on  the  basis  of  energy  balance  considerations.  For  such 
energy  analyses  it  is  convenient  to  visualize  coupled  modes  as 
shown  in  Fig.  2.  There  the  modes  are  represented  as  storage 
elements  (containing  kinetic  energies  <T.>  and  <T2>),  connected 
by  power  flow  Pi2,  supplied  with  power  inputs  and  A2  from 
sources  external  to  the  system,  and  experiencing  power  losses 
(dissipations)  D-j^  and  D2,  respectively. 
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POWER  FLOW  BETWEEN  TWO  SETS  OF  MODES 

As  has  been  discussed  earlier,  one  may  analyze  the 
vibrations  of  an  elastic  system  in  terms  of  the  vibrations 
of  the  set  of  the  modes  of  the  system.  The  interaction 
between  two  elastic  systems  may  correspondingly  be  studied 
by  investigating  the  interaction  between  two  sets  of  modes 
which  represent  the  two  systems. 

The  modes  of  an  ideal  system  vibrate  independently  of 
each  other,  as  has  been  pointed  out  previously,  and  the 
modes  of  lightly  damped  realistic  systems  tend  to  behave 
similarly.  Hence,  the  present  section  deals  with  the  power 
flow  between  two  sets  of  modes,  where  the  members  of  each 
set  are  not  coupled  to  each  other,  but  are  coupled  to  all 
members  of  the  other  set. 

The  problem  discussed  in  the  present  section  may  be 
visualized  with  the  aid  of  Fig.  3,  which  schematically  shows 
two  sets  of  modes,  one  designated  by  "a  ,  the  other  by  P  . 

It  is  assumed  that  no  power  flow  occurs  between  modes  which 
are  members  of  the  same  set7,  but  that  each  mode  of  the  a 
set  is  coupled  to  each  mode  of  the  f3  set.  It  is  desired  to 
determine  the  flow  of  energy  from  the  a  set  to  the  P  set  in 
terms  of  the  previously  studied  mode-to-mode  flows,  under 
conditions  where  each  of  the  modes  of  each  set  may  dissipate 
energy  and  may  be  supplied  with  energy  from  outside  the  system 
(as  indicated  in  Fig.  3). 

In  particular,  it  is  desired  to  determine  under  what 
conditions  one  may  obtain  relations  for  the  total  net  power 
flow  between  two  sets  of  modes  in  the  form 


p  =  <t>  T  -  T^ 

rap  ctp  a  pa  P 


(61) 


in  analogy  to  Eq.  (54),  which  pertains  to  power  flow  between 
two  single  modes.  Here  T„  and  Tp  denote  the  set-averages  of 
the  time-average  modal  kinetic  energies,  which  are  discussed 
in  more  detail  in  the  subsequent  paragraphs.  Moreover,  one 
desires  to  establish  for  coupled  sets  a  power  flow  diagram 
like  Fig.  4,  in  analogy  to  the  power  flow  diagram  of  Fig.  2, 
which  pertains  to  two  coupled  modes. 


Throughout  the  following  discussion  it  is  assumed  that 
the  power  flow  between  each  pair  of  coupled  modes  is  propor¬ 
tional  to  the  difference  between  the  time-average  modal  kinetic 
energies,  as  indicated  by  the  second  of  Eqs.  (57).  It  is  con¬ 
venient  for  the  present  purposes  to  rewirte  the  aforementioned 
expression  as 


7No  power  flow  between  these  modes  occurs  either  if  these  modes 
are  not  coupled  to  each  other,  or  if  all  of  these  modes  contain 
equal  kinetic  energies.  See  Eq.  (57). 
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(62) 


Paj,pk  "  0aj,Pk^Taj  “ 


,th 


Here  P„  .  o,  denotes  the  power  flow  from  the  j"'*1  mode  of  the 
a  set  &J*tfie  kth  mode  of  the  P  set,  and  Ta.  and  T-  denote  the 
time-average  kinetic  energies  of  these  two  modes  \aue  to 
effectively  white-noise  excitation),  in  agreement  with  Eq.  157). 
In  Eq.  (62)  and  subsequently  the  brackets  prevlously_used_to 
designate  time-averages  are  omitted  for  the  sake  of  economy  of 
notation;  however,  all  kinetic  energy  and  power  flow  quantities 
appearing  here  and  henceforth  are  intended  to  represent  time- 
average  values. 

Power  Flow  Between  Mode  Sets  in  Terms  of  Set-Average  Energies 

The  total  power  flow  Pae  from  mode  set  a  to  mode  set  P 
is  the  sum  of  all  the  individual  mode-to-mode  power  flows. 

If  set  oi  contains  Na  modes  and  set  P  contains  Np  modes,  then 
one  may  express  the  set-to-set  power  flow  as 


Na  % 

Na  Np 

Pap  =  X  X  Potj,pk  = 
j=l  k=l 

I  1 

0=1  k=l 

in  view  of  Eq.  (62) . 

The  set-average  modal 
two  sets  are  defined  as 

kinetic 

Na 

EH3 

in 

Eh3 

Te  E 

raj  "  Tpk) 


(63) 


N, 


T 


'IQ  l—> 

^  k=l 


Pk 


(64) 


If  Eqs.  (6l)  and  (63)  are  to  agree,  then  the  relation 

N~  V  V  £ 


Na  N« 

I  i 


pk^Taj  "  TPk 


J=1  k=l 


*  ~  Na  L  Taj  Np  L 

a  j-1  P  k=l 


(65) 
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must  hold  for  all  values  of  Ta,  and  To.  .  Consequently  the 
set-to-set  coupling  factors  must  obejr 


Ncc  ^ 

[Na  L  Taj]  + 

j=l  k=l 


(66) 


j=l  k=l 


0aj,Pk 


k=l 


and  the  equivalent  relations 


The  latter  expressions  have  an  interesting  interpretation. 
The  term  within  the  large  square  brackets  of  the  first  of  the 
foregoing  equations  is  the  weighted  average  of  the  0  .  ft 
coefficients  (for  a  given  value  of  k),  with  the  modal'-’,PK 
energies  Ta.  acting  as  weighting  factors.  The  term  within 
the  curved  Jbrackets  may  be  seen  to  be  the  average  (over  all 
values  of  k)  of  the  aforementioned  weighted  averages.  Thus, 

^aft  is  NqNR  tlmes  the  average  value  of  0  6  (in  the  pre¬ 

viously  indicated  sense);  this  implies  that’tfie  power  flow 
between  mode  sets  is  equal  to  the  average  mode-to-mode  flow 
times  the  number  of  mode  pairs  between  which  power  flow  occurs. 
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Modal  Energy- Independent  Coupling  Factors 


If  the  set-to-set  coupling  factors  0ao  and  0^  depend 
on  the  modal  energies  Ta.  and  To,  as  in  Eqs.  (66)*and  (67), 
then  use  of  these  factors  in  Eq.  (6l)  does  not  simplify  the 
calculation  of  power  flow  between  mode  sets.  However,  in 
some  cases  the  coupling  factors  do  turn  out  to  be  independent 
of  the  energies  Ta,  and  To  's,  and  then  Eq.  (6l)  does  present 
a  convenient  means Jfor  calculating  the  set-to-set  power  flows. 
Conditions  under  which  the  coupling  factors  tend  to  be  inde¬ 
pendent  of  the  energies  are  explored  in  the  present  section. 


Uniform  Coupling 

For  the  case  where  all  of  the  mode-to-mode  coupling 
factors  are  equal,  that  is,  where  <t>a.  =  0Q  for  all  J  and  k, 

one  finds  that  the  set-to-set  coupling  coefficients  of  Eqs.  (66) 
or  (67)  are  independent  of  the  modal  kinetic  energies  and  are 
given  by 


^ap  “  ^pa 


•Wp  *0 


(68) 


One  may  observe,  as  also  previously  mentioned,  that  here  the 
set-to-set  coupling  coefficient  is  NaNg  times  the  mode-to-mode 
coefficient,  and  that  N~No  is  the  number  of  mode-to-mode  paths 
along  which  power  can  flow  from  members  of  one  set  to  members 
of  the  other. 

Equation  (57)  indicates  that  the  assumption  of  uniform 
coupling  may  be  useful  under  some  often-encountered  practical 
circumstances.  If  all  modes  of  set  "a"  are  similar  to  each 
other  (i.e.,  have  nearly  the  same  resonance  freouency  and  the 
same  amount  of  damping),  if  all  modes  of  set  "P'r  are  similar 
to  each  other,  and  if  all  “  modes  are  similarly  coupled  to  all 
P  modes  (i.e.,  if  all  mode-to-mode  coupling  is  characterized  by 
the  same  kinds  and  magnitudes  of  coupling  elements),  then  all 
the  mode-to-mode  coupling  factors  will  be  the  same.  Thus, 

Eq.  (68)  applies  to  mode  sets  that  satisfy  the  foregoing 
similarity  criteria.  Indeed,  an  essential  step  in  the  prac¬ 
tical  application  of  the  statistical  energy  approach  consists 
of  organizing  the  set  of  all  modes  of  a  system  into  subsets 
within  which  all  modes  have  similar  characteristics,  so  that 
relations  like  Eq.  (68)  can  be  applied  to  these  subsets. 
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Uniform  Modal  Energies 

If  all  of  the  modal  kinetic  energies  in  set  "a"  are 
equal,  and  if  the  same  is  true  of  set  "P",  so  that  one  may 
write  Ta.=T  for  all  j  and  T 
or  (67)  are  round  to  reduce 


|^=Tp0  for  all  k,  then  Eqs.  (66) 


<P, 


aj,Pk 


(69) 


This  expression  indicates  that  here  the  set-to-set  coupling 
factor  is  equal  to  the  sum  of  the  mode-to-mode  factors  for 
all  coupled  mode  pairs.  (This  sum  clearly  is  equal  to  N^Ng 
times  the  average  value  of  the  mode-to-mode  coupling  factor 
between  all  "a"  and  "P"  modes.) 


In  order  to  investigate  the  circumstances  under  which 
the  assumption  of  uniform  modal  energies  is  a  realistic  one, 
one  needs  to  study  the  general  modal  kinetic  energy  To,  . 

An  expression  for  this  energy  may  be  obtained  from  a  p 
steady-state  time-average  energy  balance  relation  for  the  k— 
mode  of  set  "P",  which  relation  may  be  written  as 


or 


Na 

APk  +  X  Paj,Pk  =  Dpk 
j=l 


Na 

Apk  +  X  0aj,pk^Taj-Tpk^  =  cpk  mj^  TPk 


j=l 


(70) 


(71) 


As  indicated  in  Pig.  3,  Ag^  represents  the  power  supplied  to 
the  mode  directly  from  external  sources.  Do,  denotes  the  power 
dissipated  by  the  mode  (which  depends  on  tne  modal  viscous 
damping  coefficient  c^  and  on  the  mean  square  velocity,  and 
thus  on  the  modal  kinetic  energy  Tg^) ,  and  Pal  g^.  denotes  the 
power  supplied  from  the  j^h  mode  or  the  "a"  setpro  the  mode 
under  consideration. 
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Solution  of  Eq.  (71)  for  Tpk  yields 


which  expresses  the  kinetic  energy  Tp^  of  the  single  mode  Pk 
in  terms  of  the  modal  and  coupling  parameters  and  of  the 
energies  Ta,  of  all  of  the  modes  of  the  "a"  set.  One  may f use 
Eq.  (72)  to J investigate  under  what  conditions  all  of  the  P 
modes  possess  approximately  equal  time-average  kinetic  energies. 

If  those  terms  of  Eq.  (72)  which  contain  the  coupling 
factors  are  negligible  compared  to  the  others,  then  this 
equation  reduces  to 


T«c  “  APk  W2  °Pk 


(73) 


For  negligible  coupling  the  power  input  may  be  shown  to  be 
given  by8 


8If  coupling  effects  are  negligible,  then  the  power  input  is 
equal  to  the  power  dissipation,  or  (omitting  all  subscripts 
since  only  one  system  is  being  considered  here) 

A  =  D  =  c<V2>  =  cSI  =  vS/m  , 
in  view  of  Eqs.  (52)  and  (56). 

One  may  verify  that  one  may  obtain  the  same  result  by  direct 
calculation  of  the  power  input: 

A  =  <F(t)  V(t)>  =  /  y(a)  <F(t-a)  F(t)>  da  . 

-00 

Since  the  force  autocorrelation  function  Rf(a)  is  related  to 
the  spectrum  S  of  a  white  noise  force  (Ref.  12)  as 

<F(t-a)  F(t)>  =  Rf(a)  =  2ir  S&  6(a)  , 

where  6(a)  is  the  delta  or  unit  impulse  function,  the  above 
reduces  to  w 

A  =  2vS  f  y(a)  6(a)  da  =  TrS[y(CH-)  +  y(0-)]  =  irS/m  , 

-00 

since  y(0+)  =  1/m,  y(O-)  =  0. 
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A3k  "  r  Spk/mpk 


(74) 


9 


where  Spk  denotes  the  spectral  density  of  the  force  which 
acts  directly  on  the  3k  mode.  One  may  then  rewrite  Eq.  (74)  as 

T3k  ~  (7r/2) (Spk/Cpk^  •  (75) 

This  result  indicates  that  in  the  presence  of  small  coupling 
the  kinetic  energies  of  all  "3"  modes  will  be  equal  if  the 
excitations  are  adjusted  so  that  the  ratio  of  modal  force 
spectral  density  to  modal  damping  is  the  same  for  all  modes. 

As  a  special  case,  if  all  modes  of  the  3  set  are  equally 
excited  and  equally  damped,  then  they  will  have  equal  kinetic 
energies . 

Some  other  observations  relevant  to  uniform  distribution 
of  energies  may  be  made  from  an  analysis  of  indirectly  excited 
modes.  If  the  modes  of  the  "3"  set  are  lightly  damped  and  not 
directly  excited  (Aoj=0)  then  the  net  power  flow  to  the  jth 
mode  from  all  the  modes  of  the  "a"  set  must  be  very  small, 
and  one  may  determine  from  Eq.  (71)  or  Eq.  (72)  that 


(76) 


Equation  (76)  shows  Tpk  here  to  be  a  weighted  average  of 
all  of  the  Taj's  (with  the  coupling  coefficients  0  ,  acting 
as  weighting  factors) .  One  notes  that  the  kinetic  eftergies 
T3kOf  all  of  the  "3'1  modes  will  be  equal  if  the  sums  appearing 
in^q.  (76)  are  independent  of  the  index  k,  or  if  all  of  the 
"a"  modes  have  the  same  kinetic  energy  Ta.f  Thus,  If  the  total 
coupling  to  all  "3"  modes  is  uniform,  and/or  if  the  energy 
distribution  in  the  "a"  set  is  uniform,  then  the  energy 
distribution  in  the  indirectly  excited  ,,3"  set  will  also  be 
uniform. 


Equation  (76)  also  implies  that  the  distribution  of 
energies  Tp^  in  the  indirectly  excited  ,,3,,  set  will  generally 
be  smoother  than  that  of  the  energies  Taj  in  the  "a"  set,  if 
the  summations  do  not  depend  very  strongly  on  the  index  k. 

In  other  words,  the  Tp^'s  will  approach  equality  to  each  other 
more  closely  than  do  the  Taj's.  provided  that  the  coupling  of 
the  "3"  modes  to  all  of  the  "a"  modes  is  about  the  same  for  all 
of  the  "3"  modes. 
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One  may  thus  visualize  the  behavior  of  a  cascade  of 
lightly  damped  mode  sets,  where  only  the  first  set  is 
directly  excited,  where  each  set  is  relatively  loosely 
coupled  to  its  neighbors,  and  where  the  mode-to-mode 
coupling  for  each  pair  of  adjacent  sets  is  more  or  less 
uniform.  (This  does  not  imply  that  the  mode-to-mode 
coupling  need  have  the  same  value  for  all  set  pairs.) 

If  the  first  set  of  such  a  cascade  is  excited  so  that  an 
arbitrary  distribution  of  kinetic  energies  is  produced  in 
it,  then  the  second  set  will  attain  an  energy  distribution 
which  is  more  uniform  than  that  in  the  first,  in  accordance 
with  the  previous  paragraph.  Similarly,  the  energy  distribu¬ 
tion  in  the  third  set  will  then  be  more  uniform  than  that  in 
the  second,  and  so  on. 

Power  Flow  in  the  Ensemble  Average 

In  practice  one  rarely  deals  with  structures  whose  every 
detail  is  precisely  known.  For  example,  manufacturing  toler¬ 
ances  produce  geometric  differences,  and  variations  in  material 
properties  result  in  deviations  from  ideal  homogeneity  and 
isotropy.  Hence,  one  may  expect  to  encounter  some  differences 
in  the  properties  and  distributions  of  modes  in  nominally 
identical  structures,  and  thus  in  the  responses  of  these 
structures  to  the  same  excitation.  In  order  to  study  the 
behavior  of  "average"  structures,  one  is  led  to  consider^ the 
power  flow  that  may  be  expected  to  occur  "on  the  average" 
between  two  elastic  systems,  where  each  of  the  two  is  a  member 
of  an  ensemble,  and  where  the  members  of  each  ensemble  can 
differ  from  each  other  to  some  degree. 

Instead  of  dealing  with  only  a  single  "a"  set  and  a 
single  "P"  set  of  modes  (corresponding  to  a  single  a  elastic 
system  and  a  single  P  elastic  system),  it  is  necessary  here 
to  deal  with  collections  or  "ensembles"  of  such  sets  (corres¬ 
ponding  to  ensembles  of  elastic  systems).  The  problem  to  be 
considered  here  consists  of  determining  the  ensemble  average 
of  the  time-average  power  flow  between  a  and  P  mode  sets,  in 
terms  of  the  ensemble  average  of  the  time-average  kinetic 
energies. 

Consider  an  "experiment",  where  one  selects  at  random 
from  the  ensembles  one  set  of  the  type  a  and  one  of  the  type  P, 
subjects  each  of  the  selected  pair  of  sets  to  a  specified 
excitation,  and  measures  the  time-average  power  flow  and 
kinetic  energies.  It  is  desired  to  express  the  average  power 
flow  measured  in  a  number  of  such  experiments  in  a  form 
analogous  to  Eq.  (6l),  as 
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y 


(77) 


"  ^a(3  "  ^pa  ^p 


where  the  bar  denotes  the  average  of  the  values  obtained  in 
these  experiments  (i.e.,  the  ensemble  average)  and  the  tilde  (~) 
is  intended  to  differentiate  the  present  coupling  coefficients 
from  those  of  Eq.  (6l),  which  pertain  to  two  fully  defined  mode 
sets;  i.e.  to  a  single  "experiment". 

One  may  similarly  modify  Eqs.  (62)-(66)  to  apply  to  ensemble 
averages  and  to  the  coupling  coefficients  as  defined  by  Eq.  (77) 
and  arrive  at  expressions  analogous  to  Eqs.  (67).  Noting  that 
one  may  interchange  the  order  of  summation  and  averaging  (which 
essentially  is  another  summation),  one  obtains 


> 


(78) 


Vp 


N, 


I 

j=l 


k=l 


0<*j,Pk 


+ 


. 


The  average  of  the  product  of  two  quantities  is  not  in 

one  cannot  in 
as  one  would  like 
,  the  aforementioned 
interchange  of  the  averaging  and  multiplication  processes  is  valid 
for  random  variables  that  are  statistically  independent.  This 
means  that  0aj,Pk  Taj  =  (<^aj,Pk)  (Taj)  >  provided  that  the  probability 

of  obtaining  any  value  of  Taj  in  a  given  "experiment"  is  unaffected 

by  the  value  of  </>aj  ^  in  that  experiment,  and  vice  versa. 


general  equal  to  the  product  of  the  averages; 
general  replace  <Pa,  o'  T  by  (,<*  j  ,Pk)  (Ta , ) ,  . 
to  do  in  order  to  simplify  Eqs. (7°).  However 


The  various  modal  kinetic  energies  which  are  produced  by 
a  given  excitation  distribution  do  depend  on  the  various  mode- 
to-mode  coupling  factors  (as  one  can  readily  determine  from  any 
energy  balance  analysis),  and  thus  cannot  be  entirely  statistic¬ 
ally  independent  of  the  coupling.  However,  it  is  likely  that  a 
given  Taj.  will  depend  only  very  weakly  on  a  single  <t>a ^ 

particularly  if  many  modes  are  present  in  the  two  sets;  then 
these  two  quantities  may  be  very  nearly  independent.  If  the 
excitation  distribution  may  also  vary  from  experiment  to  experi¬ 
ment  (e.g.,  in  the  relatively  realistic  case  where  the  modal 
excitations  depend  on  the  sample  structure  selected  for  an 
"experiment")  then  statistical  independence  is  even  more  likely. 
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For  statistically  independent  mode-to-mode  coupling 
coefficients  and  modal  kinetic  energies  one  may  replace 

0aj,p}c  Taj  *>y  (0aj,|3k)  *  Taj  in  the  flrst  of  Ec*s-  ^  311(3 
perform  a  similar  replacement  in  the  second  of  these  equations. 
If  Tai  is  independent  of  the  mode  number  j  (that  is,  if  all 
modes*1 2 * * 5 of  set  a  have  the  same  kinetic  energy  in  the  ensemble 
average)  and  if  a  similar  statement  holds  about  the  P  modes, 
then  one  obtains  from  Eqs.  (78)  the  simple  result 


$ 


k=l 


aj,Pk 


(79) 


Equation  (79) >  which  is  analogous  to  Eq.  (69),  indicates 
that  the  coupling  coefficients  to  be  used  in  Eq.  (77)  are  equal 
to  the  sum  of  the  ensemble  average  (or  the  ensemble  average  of 
the  sum)  of  the  mode-to-mode  coupling  factors.  Of  course, 

Eq.  (79)  applies  only  if  the  conditions  used  in  its  derivation 
are  met,  that  is,  if  the  modal  kinetic  energies  and  mode-to-mode 
coupling  factors  are  statistically  independent,  and  if  the 
ensemble  average  modal  kinetic  energies  are  independent  of  mode 
number  (constant  within  a  type  of  set). 

Summary  and  Conclusions 

In  this  section  it  has  been  shown  that  the  steady-state 
time-average  power  flow  from  one  set  of  modes  to  another  is 
proportional  (or  approximately  proportional)  to  the  difference 
between  the  set-average  modal  kinetic  energies  of  the  two  sets, 
provided  that  both 

1)  the  mode-to-mode  coupling  satisfies  the  conditions 
discussed  under  the  heading  of  "Power  Flow  Between 
Two  Modes",  so  that  the  power  flows  between  modes 
of  the  two  sets  may  be  taken  as  proportional  to  the 
modal  kinetic  energy  differences; 

and 

2)  either  a)  the  mode-to-mode  coupling  is  (at  least 

approximately)  the  same  for  all  mode  pairs; 

or  b)  all  modes  in  a  set  have  (at  least  approxi¬ 

mately)  equal  time-average  kinetic  energies. 
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Then  the  factor  of  proportionality;  i.e.,  the  coupling 
factor,  applicable  to  set-to-set  power  flow  is  equal  to  the 
sum  of  all  the  mode-to-mode  coupling  factors  (or  to  the 
product  of  the  number  of  mode  pairs  and  of  the  average  value 
of  all  of  the  individual  mode-to-mode  coupling  factors). 

It  has  also  been  shown  that  the  ensemble  average  of  the 
steady-state  time  average  power  flows  from  one  set  of  an 
ensemble  to  a  set  of  another  ensemble  is  proportional  to  the 
difference  between  the  ensemble  average  of  the  set  average 
modal  kinetic  energies  of  the  two  sets,  -  provided  (l)  that 
the  first  of  the  previously  listed  conditions  is  satisfied, 
(2)  that  the  mode-to-mode  coupling  coefficients  and  modal 
kinetic  energies  are  statistically  independent,  and  (3)  that 
the  modal  kinetic  energies  in  the  ensemble  average  are  uni¬ 
formly  distributed  within  a  set.  The  applicable  coupling 
factor  then  is  the  ensemble  average  of  the  Siam  (or  the  sum 
of  the  ensemble  average)  of  all  the  individual  mode-to-mode 
coupling  factors. 

The  mode-to-mode  coupling  is  essentially  the  same  for 
all  mode  pairs  if  all  modes  of  the  a  set  are  similar  to  each 
other,  if  all  modes  of  the  P  set  are  similar  to  each  other, 
and  if  the  same  (or  similar)  physical  mechanisms  (e.g.,  the 
same  type  of  sound-to-structure  coupling)  couples  all  modes 
of  the  a  set  to  all  modes  of  the  (3  set. 

|(  AH  modes  in  a  set  have  equal  kinetic  energies  (i.e., 
equipartition  of  energy"  occurs  among  the  modes  of  a  set). 


either  a)  the  modes  of  the  set  are  uniformly  excited 
and  either  weakly  or  uniformly  coupled  to 
modes  of  other  sets; 

or  b)  the  modes  are  lightly  damped,  not  excited 

directly,  and  uniformly  coupled  to  a  set  of 
modes  within  which  all  modes  have  nearly  the 
same  energy. 

1x1  the  special  cases  discussed  in  this  section,  and  only 
in  these  cases,  the  average  response  and  dynamic  interaction 
of  mode  sets  can  be  reduced  to  simple  terms,  represented 
schematically  in  Pig.  4.  The  most  important  expressions 
pertaining  to  these  special  cases  are  also  summarized  in 
this  figure.  In  analyzing  the  total  interaction  of  two 
systems  by  means  of  the  statistical  energy  approach,  one 
must  generally  classify  the  modes  of  the  two  systems  into 
sets,  such  that  the  modes  within  each  set  are  sufficiently 
similar  so  as  to  satisfy  the  conditions  which  permit  this 
simple  treatment  of  set  interaction  to  be  applied. 
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LOSS  FACTORS,  RESISTANCES,  AND  MODAL  DENSITIES 

The  foregoing  development  of  the  concepts  of  the 
statistical  energy  method  was  presented  in  terms  of  power 
flow  coupling  coefficients  and  in  terms  of  the  number  of 
modes  in  a  set.  These  parameters  are  relatively  unfamiliar, 
and  generally  more  difficult  to  estimate,  calculate,  or 
measure  than  such  more  commonly  used  quantities  as  loss  factors, 
coupling  resistances,  and  modal  densities.  It  is  the  purpose 
of  the  present  section  to  present  the  relations  between  these 
more  common  parameters  and  those  used  in  the  previous  sections, 
so  that  one  may  restate  the  previously  derived  results  in  terms 
of  the  more  familiar  quantities. 

The  relations  between  these  two  sets  of  parameters  are 
simple,  so  that  restatement  of  the  previous  results  in  terms 
of  the  more  common  quantities  involves  only  trivial  algebraic 
manipulations.  Hence,  no  restated  results  are  given  here. 

Some  applications  of  the  relations  given  in  the  present  section 
are  illustrated  in  the  subsequent  one. 


Loss  Factors 

The  loss  factor  is  classically  defined  on  the  basis  of  a 
single-degree-of-freedom  system  oscillating  sinusoidally  at 
resonance  in  the  steady  state.  The  definition  of  the  loss 
factor  t)  may  be  stated  as 

t)  =  Ed/2rrE  ,  (80) 

where  E.  denotes  the  energy  dissipated  per  cycle  by  the  system 
and  E  denotes  the  "energy  of  vibration",  which  is  equal  essen¬ 
tially  to  the  total  (kinetic  plus  potential)  energy  at  any  time 
(and  is  nearly  constant),  to  the  time-wise  maximum  kinetic  energy, 
and  also  to  the  time-wise  maximum  potential  energy  of  the  system 
(Refs.  13,  14). 

One  may  express  the  loss  factor  in  terms  of  the  time-average 
power  D  dissipated  by  the  system,  since  D  =  E,  o> /2ir,  where  co 
denotes  the  radian  natural  frequency  of  the  system  so  that  o> 72* 
represents  the  natural  frequency  in  cycles/time.  One  obtains 


D  _  D 

4  "  5THE  "  25TT 

o  o 


(81) 
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where  T  denotes  the  time-average  kinetic  energy  and  T=E/2 
holds  for  the  postulated  oscillations.  Since  for  a  simple 
mass-spring-dashpot  system  D  =  c<v^>  and  T  =  <v2>m/2,  one 
may  readily  verify  that  for  such  a  system 


1  = 


c  _  c 
^o  Jtol 


(82) 


where  c  denotes  the  critical  viscous  damping  coefficient  of 
the  system.  Equation  (82)  involves  the  viscous  damping 
coefficient  c,  and  hence  applies  only  for  damping  of  the 
viscous  type;  Eq.  (82)  thus  is  less  general  than  Eqs.  (80) 
and  (8l),  which  are  not  restricted  to  any  particular  type  of 
energy  dissipation  mechanism.  Since  the  loss  factor  applies 
for  all  damping  mechanisms,  it  is  generally  preferable  to 
more  restricted  measures  of  damping. 

A  broadened  interpretation  of  Eq.  (8l)  is  often  used  to 
define  (and  measure)  loss  factors.  In  this  interpretation 
a)  is  taken  as  the  center  frequency  of  a  band  encompassing 
tne  system  natural  frequency,  E  is  taken  as  the  system  energy 
in  that  band,  and  D  is  taken  as  the  power  dissipated  by  or 
(transferred  away  from)  the  system  in  the  same  band.  This 
interpretation  is  generally  not  needed  for  dealing  with  a 
single-degree-of -freedom  system,  but  is  useful  when  one  is 
concerned  with  a  collection  of  such  systems,  -  i.e.,  with  a 
set  of  modes.  For  a  set  of  modes,  D  of  Eq.  (8l)  is  taken  as 
the  total  power  dissipated  by  all  modes  and  T  as  total  kinetic 
energy  of  all  modes.  For  a  set  a  containing  N  individual 
modes,  the  total  power  dissipation  D_  and  total  time-average 
kinetic  energy  TTa  are  given  by  1 


J=1 
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where  q  .  denotes  the  loss  factor  of  the  aj  mode,  as  defined 
by  Eq.  TJ8l)  and  Ta  denotes  the  set-average  kinetic  energy, 
as  in  Eq.  (64).  The  loss  factor  qa  of  set  a  then  obeys 


D, 


2co  T, 


'Ta 


Ta 


2cjd  N  T 
6  aAa 


(84) 


where  now  co  is  defined  as  the  center  frequency  of  a  band 
which  encompasses  all  modal  resonances.  The  last  form  of 
the  above  equation  indicates  that  the  set  loss  factor  qa  is 
a  weighted  average  of  the  individual  modal  loss  factor  qa,, 
with  the  modal  kinetic  energies  TaJ  acting  as  weighting 
factors. 


Coupling  Loss  Factors 

Mechanical  energy  may  be  lost  from  a  given  system  not 
only  by  being  dissipated  (i.e.,  being  changed  to  thermal 
energy),  but  also  by  being  transferred  to  another  system. 
Energy  conducted  from  a  system  is  lost  to  it,  just  like 
energy  dissipated  by  it.  Hence,  it  is  useful  to  define  a 
loss  factor  for  energy  conduction  analogous  to  that  for 
energy  dissipation. 

It  is  convenient  to  define  a  coupling  loss  factor  qa,  ok 
of  mode  a.1  to  describe  the  time-average  power  flow  Pa  .  pk 
from  mode  aj  to  mode  Pk  in  analogy  to  Eq.  (8l)  by 


^j,Pk 


Fo.1tPk 


2  0) 

o 


(85) 


The  condition  To.  =0  is  imposed  in  order  to  make  this  loss 
factor  independent  of  the  modal  kinetic  energies.  In  view 
of  Eq.  (62),  one  finds  that  the  coupling  loss  factor  is 
related  to  the  mode-to-mode  coupling  factor  as 


naj,Pk 


(86) 
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Eq.  (86)  indicates  that  the  same  coupling  loss  factor  applies 
for  power  flow  in  either  direction  (that  is,  q 
for  those  various  previously  discussed  cases  lA1 
the  same  coupling  factor  applies  for  power  flow  in  either 
direction  <0aJ;Pk  =  *pk>ctJ). 


One  may  also  define  an  analogous  loss  factor  to  account 
for  power  flow  from  one  set  of  modes  (a)  to  another  (p).  The 
coupling  loss  factor  T)ap  of  set  a  is  defined  as 


•ap 


rap 


2u)o  V 


v° 


(87) 


where  a>  denotes  the  center  frequency  of  a  band  which  encompasses 
all  modal  resonances  of  both  sets.  In  view  of  Eqs.  (6l)  and  (83), 
one  may  rewrite  the  foregoing  as 


^ap 

Vp  -  ga>0  Na  Ta 


ap 


2  %  Na 


(88) 


Here  it  is  important  to  note  that  because  of  the  appearance  of 
Nqt  in  Eq.  (88)  the  coupling  loss  factor  q6a  of  set  P  (for  power 
flow  from  P  to  a)  is  generally  not  the  same  as  the  loss  factor 
Tlap  of  set  a  (for  Power  flow  from  a  to  P),  even  if  the  coupling 
coefficients  0ao  =  0^  are  the  same  for  both  flow  directions 
(as  /oo\he  PrevionsiF’ analyzed  cases).  Instead,  in  view  of 
Eq.  (88)  and  an  expression  analogous  to  it. 


^ap  Na  ”  ^pa  Np  “  ^ap^2  % 


(89) 


Loss  and  Coupling  Resistances 

In  analogy  to  electrical  circuits,  in  which  the  time- 
average  power  dissipation  is  given  by  <I2>R  ,  where  I  denotes 
current  and  R0  electrical  resistance,  one  mly  define  a  mechanical 
or  acoustic  resistance  R  of  an  elastic  system  so  that  the  power 
loss  D  from  the  system  is  given  by 

D  =  R  <V2>  .  (90) 
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For  a  viscously  damped  single-degree-of-freedom  system  the 
resistance  is  equal  to  the  viscous  damping  coefficient  c, 
of  course.  In  view  of  Eq.  (82)  one  may  relate  this  resistance 
to  the  loss  factor  q  and  obtain 

R  =  c=  qma>o  (91) 

for  a  system  of  mass  m  and  natural  frequency  o)q. 

One  may  define  a  coupling  resistance  1^,  which  describes 
the  power  loss  of  a  mode  aj  to  a  mode  Pk  anaidgously  to  Eq.  (90) 
by 


<V?,>  = 


(92) 


In  view  of  Eqs .  (85)  and  (30)  one  may  determine  that  the  coupling 
resistance  is  related  to  the  corresponding  loss  factor  as 

=  "kj  “o  *  ^ 


Since  this  expression  is  seen  to  involve  the  modal  mass  n^., 
one  observes  that  in  general  although  theJ 

corresponding  two  loss  factors  mdy oe  equdlr 

The  loss  resistance  applicable  to  a  set  a  of  modes 
should  describe  the  total  power  loss  according  to  an  equation 
like  Eq.  (90).  The  choice  of  what  one  defines  as  the  mean 
square  velocity  is  somewhat  arbitrary;  it  turns  out  that  in 
dealing  with  a  set  of  modes  of  an  elastic  system  of  mass  m^ 
it  is  convenient  to  choose  <V£>  =  2  T^/n^.  Then  one  obtains 


*a  =  '2  T^  =  %  %  *  (94) 

in  view  of  Eq.  (84).  The  above  expression  is  analogous  to 
Eqs.  (91)  and  (93),  except  that  here  represents  a  system 
mass,  rather  than  a  modal  mass.8 9 


8With  the  mode  shape  normalization  indicated  in  Eq.  (7),  the 

modal  mass  i£  the  same  as  the  system  mass.  But  other  normal¬ 
izations  may  be  used  legitimately,  and  with  those  this  equality 
would  not  hold. 
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By  analogy  to  Eq.  (94)  one  may  define  a  coupling 
resistance  R^  to  account  for  the  loss  of  power  from  a  set 
of  modes  o.  topa  set  of  modes  P  as 


_  Pap  ma 
2 


(95) 


and  determine  by  use  of  Eqs.  (87)  and  (88)  that 


0ctp  "Vx 

'xxp  “  ^ap  "Vx  %  “  2  Na 


(96) 


Again,  R^  and  R^a  are  generally  different,  even  If  0ap=0pa. 
Modal  Densities 

One  usually  studies  the  interactions  of  elastic  systems 
on  a  band-by-band  basis,  by  determining  the  power  flow  from 
modes  of  one  system  which  have  their  resonances  in  a  given 
frequency  band  to  a  set  of  modes  of  another  system  which  have 
their  resonances  in  the  same  band.  One  may  in  general  choose 
any  convenient  bandwidth,  and  one  may  often  wish  to  use 
different  bandwidths  for  different  calculations.  The  number 
of  modes  which  fall  within  the  various  sets  depends  on  the 
bandwidths  one  chooses,  so  that  it  is  useful  to  introduce 
the  concept  of  modal  density,  in  terms  of  which  the  dependence 
of  some  of  the  previous  results  on  bandwidth  can  be  stated 
more  explicitly. 

In  addition,  it  turns  out  that  modal  densities  for  many 
systems  can  be  estimated  readily;  thus,  use  of  modal  densities 
increases  the  practical  utility  of  the  analysis  method  dis¬ 
cussed  in  this  report. 

The  modal  density  of  an  elastic  system  is  defined  as  the 
average  number  of  modes  per  unit  frequency  interval.  If  a 
system  exhibits  Na  modes  whose  resonances  fall  within  a  fre¬ 
quency  interval  Aoo,  then  the  modal  density  of  the  system  at 
the  center  frequency  00  of  the  interval  is  defined  as 


na  (a))  =  Na/Aa> 


(97) 
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For  most  purposes  it  suffices  to  consider  na(<o)  as  given  by 
Eq.  (97)  to  be  a  "smooth"  function  of  frequency.  Such  smooth¬ 
ness  may  be  attained  by  restricting  the  applicability  of 
Eq.  (97)  to  cases  where  the  frequency  interval  of  interest 
encompasses  several  modes.10 

The  modal  densities  of  systems  may  be  determined  from 
their  "frequency  equations",  i.e.,  from  the  equations  which 
give  the  system  resonance  frequencies  as  a  function  of  the 
system  parameters.  For  simple  systems  these  calculations  may 
be  carried  out  without  great  difficulty,  as  illustrated  in 
Ref.  5. 


Table  I  lists  the  modal  densities  of  some  uniform  elastic 
systems,  as  obtained  from  Ref.  5  and  from  calculations  based 
on  Ref.  2.  The  expressions  listed  in  the  table  apply  strictly 
only  for  frequencies  which  are  considerably  above  the  system 
fundamental  (by  perhaps  at  least  two  octaves),  where  boundary 
conditions  have  no  important  effects;  however,  these  expressions 
generally  also  provide  reasonable  estimates  for  the  modal  densi¬ 
ties  at  frequencies  only  slightly  above  the  fundamental. 

The  modal  density  of  a  composite  system  is  approximately 
equal  to  the  stun  of  the  modal  densities  of  the  component 
systems.  Thus,  for  example,  the  modal  density  of  a  plate  with 
attached  beams  is  roughly  equal  to  the  modal  density  of  the 
plate  by  itself,  plus  the  modal  densities  of  all  of  the  beams 
by  themselves.  This  additive  property  of  modal  densities  has 
not  been  validated  analytically,  but  appears  plausible  on  the 
basis  of  some  limited  experimental  evidence  and  the  following 
reasoning  (Ref.  17). 


lo0therwise,  small  changes  in  ou  might  result  in  large  stepwise 
changes  in  na.  For  example,  if  A<n  is  small,  it  might  encom¬ 
pass  one  or  more  resonances  for  one  center  frequency,  but  none 
for  another  frequency  near  the  previous  one.  The  assumption 
of  smoothness  generally  introduces  no  important  errors  in 
estimates  of  average  responses  in  broad  frequency  bands. 

(The  errors  due  to  this  assumption  usually  are  much  smaller 
than  those  introduced  by  other  approximations  and  estimates 
one  usually  must  make  in  order  to  apply  the  statistical 
energy  approach  to  practical  problems.)  However,  one  may 
need  to  concern  oneself  with  the  fluctuations  in  na  in  dealing 
with  more  advanced  problems  concerning  the  deviations  of 
response  from  the  broad-band  average. 
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MODAL  DENSITIES  OP  SOME  UNIFORM  SYSTEMS* 
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*See  next  page  for  definitions  of  symbols. 


Symbol  Definitions  for  Table  I 


A  cross-section  area 

A„  surface  area 
cQ  acoustic  wave  velocity 

C£  longitudinal  wave  velocity 

cm  membrane  wave  velocity 

cs  string  wave  velocity 

cT  torsional  wave  velocity 

D  plate  rigidity 

E  Young ' s  modulus 

G  shear  modulus 

h  thickness 

I  centroidal  moment  of  inertia  of  A 

J  polar  moment  of  inertia  of  A 

K  torsional  constant  of  A 

L  length 

S  membrane  tension  force/unit  edge  length 

T  string  tension  force 

VQ  volume 

radius  of  gyration  of  A 

Kp  radius  of  gyration  of  plate  cross  section 

v  Poisson's  ratio 

to  frequency  (radians/time) 

p  material  density 
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The  number  of  modes  in  a  lumped  parameter  system  (e.g., 
one  consisting  of  ideal  springs,  masses,  and  dasnpots)  is 
equal  to  the  number  of  degrees  of  freedom.  If  one  inter¬ 
connects  two  such  systems,  then  the  total  number  of  degrees 
of  freedom  of  the  composite  system  is  equal  (or  very  nearly 
equal,  depending  on  the  manner  of  interconnection)  to  the 
sum  of  the  number  of  degrees  of  freedom  of  the  individual 
systems.  Thus,  the  number  of  modes  of  the  combined  system 
is  equal  to  the  sum  of  the  numbers  of  modes  of  the  two  con- 
stitutent  systems. 

However,  some  additional  considerations  are  required, 
since  the  modal  resonances  of  the  combined  system  generally 
occur  at  different  frequencies  than  the  resonances  of  the 
individual  systems.  If  one  interconnects  two  single-degree- 
of-freedom  systems  having  different  resonant  frequencies, 
then  one  obtains  a  two-degree-of-freedom  system  which  has  one 
resonance  below  the  lower  of  the  two  individual  system  reso¬ 
nances,  and  one  which  is  above  the  higher  of  the  individual 
resonances.  In  other  words,  the  interconnection  in  effect 
shifts  the  resonances  both  upward  and  downward.  If  two  multi 
modal  systems  are  interconnected,  it  is  likely  that  on  the 
average  as  many  resonances  are  shifted  into  a  fixed  frequency 
band  as  are  shifted  out  of  it,  so  that  the  total  number  of 
modes  in  the  band  is  unchanged  by  the  interconnection  process 
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SOME  SPECIAL  RESULTS  FOR  TWO  COUPLED  MODE  SETS 

For  two  coupled  mode  sets,  where  one  of  the  sets  is  not 
directly  excited,  use  of  the  statistical  energy  method  is 
particularly  fruitful  and  can  lead  to  valuable  qualitative 
insights,  as  well  as  to  quantitative  results.  The  present 
section  summarizes  general  results  applicable  to  the  afore¬ 
mentioned  case,  points  out  some  conclusions  one  may  draw 
from  them,  and  illustrates  application  of  these  results  to 
some  special  cases. 


Vibrations  of  Indirectly  Excited  Mode  Sets 

Consider  two  coupled  sets  of  modes,  as  represented 
schematically  in  Fig.  4,  and  assume  that  Ao=0,  -  that  is, 
that  no  energy  is  supplied  to  mode  set  P  except  via  Its 
coupling  to  mode  set  a.  In  the  steady  state  the  total 
power  Dtio  dissipated  by  set  p  must  be  equal  to  the  power  Pa6 
supplieato  ifc  from  set  a.  If  the  modes  of  each  set  are  p 
sufficiently  "similar",  in  the  sense  indicated  in  the  previous 
chapter,  so  that  one  may  validly  use  Eq.  (6l)  with  ^gp^p a  to 
describe  Pa0,  then  the  energy  balance  for  set  p  requires  that 

=  2  VWe  ’  (98) 

in  view  of  Eq.  (84).  If  one  solves  this  expression  for  To  one 
obtains  a  relation  between  the  average  modal  kinetic  energies 
of  the  two  sets  which  may  be  expressed  as 


(99) 


where  Eq.  (89)  has  been  used  to  replace  by  the  coupling 
loss  factor11  tj^. 


xlIt  should  be  noted  that  Eq.  (99)  involves  the  loss  factor 

which  pertains  to  energy  flow  from  P  to  whereas  the  actual 
energy  flow  must  be  in  the  other  direction  if  set  P  is  not 
directly  excited. 
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As  indicated  in  the  discussion  pertaining  to  Eqs.  (17) 
and  (18),  the  space-  and  time-average  mean  square  velocity 

<Va>  associated  with  Nq  modes  of  an  elastic  system  of  mass 

ma  is  related  to  the  total  kinetic  energy  T  of  these  modes 
as  ia 

TTa  =  NaTa  =  k  ma<^a>  *  (100) 

If  one  uses  this  expression  and  an  analogous  one  for  set  P 
to  eliminate  the  average  modal  energies  T  and  TR  from 
Eq.  (99)  >  one  obtains  the  result 


=  .  V  _  Vm>  .  V 

<V®>  Na  me  W'p  mp  '  V+,)|3 


The  last  form  of  this  expression  is  obtained  by  use  of  Eq.  (97) 
and  is  based  on  the  assumption  that  the  N  and  NR  modes  correspond 
to  the  same  frequency  band  Aa>.  up 

Equation  (101)  is  a  surprisingly  simple  and  useful  result, 
which  permits  one  to  calculate  the  mean  square  velocity  of  the 
indirectly  excited  elastic  system  from  that  of  the  directly 
excited  one  if  one  knows  the  masses  of  the  two  systems,  their 
modal  densities  (e.g.,  from  Table  I)  in  the  various  frequency 
bands  of  interest,  and  the  coupling  and  dissipation  loss  factors. 

A  further  simplification  occurs  for  those  cases  where  qa«qft  , 
since  then  Eq.  (101)  reduces  to  a  relation  which  does  not  involve 
the  loss  factors.  Applications  of  Eq.  (101)  to  a  beam  coupled  to 
a  plate  and  to  two  coupled  plates,  and  related  experimental 
results  showing  generally  good  agreement  with  theoretical  pre¬ 
dictions,  are  presented  in  Ref.  18. 

The  grouping  of  the  loss  factor  terms  appearing  in  Eqs.  (99) 
and  (101)  permits  one  to  make  some  observations  concerning  the 
effectiveness  of  adding  damping  to  indirectly  excited  systems. 

If  the  damping  of  the  p  system  is  small  initially,  that  is,  if 
^p<<T^PaJ  then  added  damping  will  reduce  the  system  vibrations 
—2 

<Vp>  only  if  this  added  damping  (or  increase  in  q^)  is  significant 

as  compared  to  qga.  This  observation  may  explain  why  the  addition 
of  some  damping  to  a  structure  (system  P)  excited  by  sound  in  an 
acoustic  space  (system  a)  may  have  little  effect  on  some  compo¬ 
nents  of  the  structural  vibrations. 
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On  the  other  hand,  Eq.  (lOl)  indicates  that  the  mean 
square  velocity  of  the  indirectly  excited  system  varies 
approximately  inversely  with  the  loss  factor  t|r  in  those 
cases  where  Thus,  added  damping  is  effective  in 

reducing  vibrations  of  indirectly  excited  systems,  provided 
that  the  resulting  dissipation  loss  factor  exceeds  the 
coupling  loss  factor. 


From  Eq.  (99)  or  (98)  one  may  conclude  that  Tg==T  if 
Hr  »t1r*  That  is,  if  the  power  dissipation  of  setpP  is 
negligible,  then  so  is  the  power  flow  from  set  a  to  set  0, 
and  the  average  kinetic  energy  Tg  of  the  indirectly  excited 
set  will  approach  that  of  the  directly  excited  set.  Thus, 
if  T| _  »t|r,  then  the  two  mode  sets  will  approach  "modal 
energy  equilibrium"  or "modal  energy  equipartition" . 


Effect  of  Coupling  on  Vibrations  of  Directly  Excited  Mode  Sets 

The  present  discussion,  like  that  of  the  previous  section, 
deals  with  two  coupled  sets  of  modes,  whose  interaction  may  be 
visualized  by  means  of  the  lower  portion  of  Fig.  3.  Here,  as 
in  the  previous  section,  no  energy  is  assumed  to  be  supplied  to 
set  0  from  outside  the  system;  i.e.  A^=0.  However,  unlike  the 
previous  section,  which  focused  on  the  indirectly  excited  set  0, 
the  present  section  is  concerned  with  the  directly  excited  set 

Conservation  of  energy  requires  that  the  time-average  power 
A  supplied  to  set  a  in  the  steady  state  be  equal  to  the  total 
power  dissipated  by  both  mode  sets.  As  previously  pointed  out 
the  power  dissipated  by  set  0  must  be  equal  to  the  power  supplied 
to  it  from  set  a.  Thus,  one  may  write 


Aa  =  2cDoT^aNaTa  +  0ap(Ta"Tp) 

-  a»oTa[v\»  +  SET  C1  -  l£)] 


(102) 


If  one  eliminates  Tft  by  substituting  Eq.  (99),  and  if  one  intro¬ 
duces  the  coupling  loss  factors  by  use  of  Eq.  (89),  tnen  one  may 
find  that  the  foregoing  results  in 
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(103) 


"a 


2%TaNa 


=  n 


'a 


=  n 


^ap 


app 


a 


*pa 


+  *1 


P 


The  "apparent  loss  factor"  Tja  Introduced  In  the  fore- 

app 

going  equation  corresponds  to  the  value  of  the  dissipation  loss 
factor  which  one  would  ascribe  to  set  a  (on  the  basis  of  measure¬ 
ments  performed  on  set  a)  if  one  were  not  aware  that  this  set  is 
coupled  to  set  P.  In  order  to  determine  the  loss  factor  qa  of 
set  a  in  absence  of  coupling,  one  would  turn  to  Eq.  (84)  and 
determine  the  values  of  the  various  parameters  Na,  Ta,  that 

appear  there.  Since  the  directly  supplied  average  power  Aa  must 
be  equal  to  the  total  dissipated  (in  absence  of  coupling  to 
other  mode  sets),  one  may  replace  D™  by  Aa;  then  the  ratio 
appearing  in  Eq.  (84)  becomes  precisely  that  which  occurs  on  the 
left-hand  side  of  Eq.  (103). 

It  is  evident  from  Eq.  (103)  that  the  apparent  loss  factor 
T]a  is  never  smaller  than  the  actual  dissipation  loss  factor  n 
app  a 

Thus,  if  one  is  unaware  that  a  system  whose  loss  factor  one  is 
measuring  is  coupled  to  another  system,  then  one  always  obtains 
a  loss  factor  value  which  is  too  large.  The  error  is  insignifi¬ 
cant,  however,  if  the  coupling  is  poor  (RaajRar,  «  0)  and/or  if 
the  coupled  set  is  relatively  lossless  (Tio«qT).  Effects  of 
energy  conduction  on  loss-factor  measurements  are  discussed  in 
some  more  detail  in  Refs.  19  and  20. 

From  Eq.  (103)  one  may  also  deduce  that  the  response 

<va>  =  2^aTa//ina  to  a  given  excitation  Aa  is  not  controlled  by 

the  actual  dissipation  loss  factor  qa  (as  in  the  uncoupled  case), 
but  rather  by  the  apparent  loss  factor  tj  .  Thus,  the  coupled 

app 

indirectly  excited  system  P  increases  the  effective  damping  of 
the  directly  excited  system  a;  -  a  conclusion  which  is  intuitively 
obvious",  since  system  P  here  serves  to  extract  energy  from 
system  a. 

Interaction  of  Sound  and  Structures 

Some  of  the  previous  results  of  this  section  are  particularly 
useful  for  estimating  the  sound  produced  in  an  acoustic  volume  due 
to  the  vibration  of  a  structure,  or  for  predicting  the  structural 
vibrations  induced  by  a  diffuse  sound  field. 
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Consider  the  case  where  a  structure  S  is  directly  excited 
by  random  forces  acting  on  it,  and  where  these  structural 
vibrations  produce  noise  in  a  room  R  in  which  the  structure 
is  located.  Here  the  structural  modes  in  a  frequency  band  of 
interest  correspond  to  a  directly  excited  set  a,  the  modes  of 
the  acoustic  space  R  correspond  to  the  indirectly  excited  set  P. 

The  total  energy  ER  in  a  diffuse  acoustic  field  within 
fluid  volume  ^in  a  frequency  band  Aoo  is  related  to  the  mean 
square  acoustic  pressure  <pO  measured  in  the  same  frequency 
band  as  (p.  95  of  Ref.  5) 

er  -  <P2>  VPaca  <10i,) 


where  p  denotes  the  density  and  c 
fluid.  aThis  total  energy  is  on  the 


the  sound  velocity  of  the 
average  half  potential  and 


half  kinetic,  hence  the  average  total  kinetic  energy  TTR  is  one 
half  of  the  total.  The  average  modal  kinetic  energy  TR  may  be 
obtained  by  dividing  Ttr  by  nR.Aoo,  the  number  of  modes  in  the 
interval.  With  the  moaSl  density  expression  given  in  Table  I(P.4c) 
for  acoustic  volumes,  one  thus  finds  the  average  kinetic  energy 
of  an  acoustic  mode  of  a  room  to  be  given  by 


TR  = 


R 

2  nR  Aa) 


^  ca  <p2> 
2  Aoj 

Pa  ® 


(105) 


If  one  applies  Eq.  (100)  to  the  structure  under  considera¬ 
tion  and  replaces  the  number  of  modes  Nq  by  n  -Au)  according  to 
Eq.  (97),  then  one  finds  that  the  average  kinetic  energy  of  a 
structural  mode  is  given  by 


TS  2  n  -Aco 


mS<a3> 

p 

2  oo  ns  Aco 


(106) 


p  _ 2  2  2 

where  <v  >  represents  the  mean  square  velocity  and  <a  >=oj  <7  > 

the  meanssquare  acceleration  of  the  structure  measured  in  the 

frequency  band  Acd  with  center  frequency  co. 

Introduction  of  Eqs.  (105)  and  (106)  into  Eq.  (99),  with 
R  taking  the  place  of  the  indirectly  excited  system  P  and  S 
that  of  the  directly  excited  system  ot,  leads  to  the  following 
relation  between  the  mean  square  values  of  the  acoustic  pressure 
and  the  structural  acceleration: 
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(107) 


<p2>  Pa  ms  ^RS 
<af>  2Tr2cans  nRS+T1R 


In  most  cases  the  room-to— structure  coupling  loss  factor  tj 
may  be  expected  to  be  much  smaller  than  the  dissipation  lo 
factor  of  the  room.  With  this  inequality,  and  with 


nRS 


^SR  ns  =  ns  RSR 
nR  nR  ms  00 


(108) 


obtained  from  Eqs.  (89),  (97),  (96),  one  finds  that  one  may 
rewrite  Eq.  (107)  as 


<P2>  =  pa  ^SR  =  fsR  #  Pa  ca 

<5s>  z»2ca  cd  nRi)R  V0cd3 


(109) 


where  the  last  form  is  obtained  by  substituting  for  n  the 
appropriate  expression  from  Table  I.  The  structure-to-room 
coupling  resistance  RqR  corresponds  to  the  "radiation 
resistance"  which  acts  on  the  structure;  i.e.  to  the  resistance 

which  governs  the  loss  of  energy  by  the  structure  due  to  acoustic 
radiation  from  it. 


One  may  also  consider  a  situation  in  which  an  acoustic 
volume  R  is  directly  excited  by  a  noise  source,  and  where  the 
resulting  acoustic  pressures  cause  an  (indirectly  excited) 
structure  to  vibrate.  Here  the  volume  excites  the  structure; 
i.e.  the  roles  of  the  two  systems  here  are  interchanged  from 
the  previously  discussed  case.  Equations  (104)-(106)  still 
apply,  but  now  the  structure  S  corresponds  to  the  indirectly 
excited  set  P  and  the  room  R  to  the  indirectly  excited  set  cx 
Here,  use  of  Eq.  (99)  yields  a  result  which  one  may  write 


<a?>  0  c  n 

— I—  =  2v2  -A  _s  ^ 

<P2>  pa  ms 


(110) 


where 
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(Ill) 


nSR  ^R 

nSR+T1s  RSR+RS 


and  the  last  form  of  Eq.  (ill)  follows  from  Eq.  (96). 

The  results  given  in  Eqs.  ( 109) -(ill)  first  appeared  in 
Ref.  10.  This  reference  also  describes  some  experimental 
results,  which  are  in  good  agreement  with  theoretical 
predictions . 

A  word  of  caution  is  required,  however,  concerning  the 
applicability  of  the  results  of  this  section.  These  results 
are  derived  from  Eqs.  (98)  and  (99)*  which  in  turn  are  based 
on  the  assumption  that  a  simple  power  flow  relation  holds  for 
the  two  mode  sets  considered.  Thus,  Eqs.  (l09)-(lll)  should 
only  be  applied  to  cases  where  all  structural  modes  in  the 
frequency  band  of  interest  are  uniformly  coupled  to  all  room 
modes  in  the  band.  If  differently  excited  or  differently 
coupled  modes  occur  in  the  band,  one  must  divide  them  into 
groups  of  like  modes,  and  apply  the  foregoing  results  to  each 
group  separately,  of  course  using  properly  modified  modal 
density  expressions,  as  needed.  Differences  in  sound-structure 
coupling  and  calculations  taking  these  into  account  are  dis¬ 
cussed  in  Ref.  5. 
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CONCLUDING  REMARKS 


It  is  hoped  that  this  report  will  introduce  the  reader 
to  the  basic  concepts  involved  in  the  statistical  energy- 
analysis  of  vibrating  systems,  so  that  he  may  put  this  approach 
to  use  with  some  confidence.  It  was  endeavored  to  describe 
the  most  important  aspects  of  this  analysis  approach  in 
considerable  detail,  in  order  to  give  the  reader  an  under¬ 
standing  of  its  range  of  validity.  Of  course,  not  all 
aspects  of  available  results  pertinent  to  this  approach 
could  be  covered  within  the  time  and  space  allotted  to  this 
report;  the  reader  interested  in  particular  topics  may  do 
well  to  consult  some  of  the  references  given  here,  as  well 
as  some  of  the  background  information  cited  in  these  refer¬ 
ences. 

Although  the  results  presented  in  this  report  are  limited 
to  linear  coupling,  a  study  of  nonlinearly  coupled  modes  has 
been  published  recently  (Ref.  21).  This  study  shows  that  the 
net  energy  flow  between  nonlinearly  coupled  modes  also  occurs 
in  the  direction  of  the  energy  gradient,  and  that  no  energy 
flows  if  the  modal  energies  are  equal. 

The  present  report  has  dealt  essentially  only  with  two 
coupled  modes,  or  with  two  coupled  mode  sets,  but  not  with 
more  sets  or  modes  in  cascade.  For  loose  coupling,  extension 
of  the  two-set  results  to  multiset-systems  can  readily  be 
accomplished.  Such  extensions  and  applications  of  them  appear 
in  Refs.  22-24,  and  a  complete  electric-circuit  analogy  for 
energy  flow  in  multiset -systems  is  developed  in  Ref.  25. 

As  yet  unpublished  calculations  indicate  that  the  power  flow 
between  any  two  adjacent  modes  of  three  linearly  coupled 
modes  in  cascade  is  proportional  to  the  modal  energy  difference, 
even  if  the  modes  are  well  coupled;  however,  here  the  factors 
of  proportionality  generally  involve  properties  of  all  three 
modes  and  thus  differ  from  the  two-coupled-mode  case. 

Calculations  on  the  basis  of  the  statistical  energy 
approach  of  variances  and  confidence  limits  in  addition  to 
the  commonly  computed  averages  of  the  dynamic  variables,  are 
presented  in  Ref.  18.  Illustrations  of  applications  of  the 
statistical  energy  approach  and  discussions  of  the  problems 
which  have  been  studied  by  means  of  it  appear  in  Refs.  17-19, 
22-28. 


At  present  the  usefulness  of  the  statistical  energy 
approach  in  predicting  vibration  responses  under  practical 
circumstances  is  limited  by  the  amount  of  information  that  is 
available  concerning  coupling  coefficients  and  loss  factors. 
Approximate  determination  of  sound-to-structure  coupling 
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coefficients  (or  resistances)  often  is  not  very  difficult, 
as  indicated  in  Ref.  5,  but  calculations  of  structure-to- 
structure  coupling  coefficients  may  involve  considerable 
effort  (e.g.,  see  Ref.  18).  A  reasonable  amount  of  infor¬ 
mation  on  the  dissipative  loss  factors  of  materials  and 
simple  structures  is  available  in  the  technical  literature 
(e.g..  Ref.  13),  but  the  loss  factors  of  built-up  structures 
can  still  not  be  estimated  with  much  confidence  -  although 
some  progress  is  being  made  (Ref.  29). 

The  statistical  energy  approach  is  a  uniquely  simple 
and  powerful  tool  for  providing  one  with  a  qualitative 
understanding  of  the  most  important  aspects  of  the  vibra¬ 
tions  of  complex  systems.  When  combined  with  other, 
analytically  or  experimentally  derived,  information 
concerning  the  energy  transport  and  dissipation  parameters 
involved  in  a  given  case,  this  approach  also  permits  one 
to  arrive  readily  at  quantitative  answers  to  complex 
vibration  problems. 
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FIG. 2  POWER  FLOW  DIAGRAM  FOR  TWO  COUPLED  MODES 
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FIG. 3  TWO  COUPLED  SETS  OF  MODES 
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See  next  page  for  definitions  of  symbols. 
See  pp.  39,  40  for  summary  of  restrictions 
under  which  the  above  relations  apply. 


FIG. 4  POWER  FLOW  BETWEEN  TWO  SETS  OF  MODES 
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SYMBOL  DEFINITIONS  FOR  FIG.  4 


The  following  abbreviations  are  used  in  this  list  of  symbols: 

t.a.  time-average  tot.  total 

s.a.  set-average  m.s.a.  mode(s)  of  set  a 


Symbols 

Definitions 

ATa 

t.a.  tot.  power  supplied  to  all  m.s.a. 

DTa 

t.a.  tot.  power  dissipated  by  all  m.s.a. 

TTa 

t.a.  tot.  kinetic  energy  of  all  m.s.a. 

Aa 

s.a.  of  t.a.  power  supplied  to  single  m.s.a. 

Da 

s.a.  of  t.a.  power  dissipated  by  single  m.s.a. 

T 

a 

s.a.  of  t.a.  kinetic  energy  of  single  m.s.a. 

N 

a 

number  of  modes  in  set  a 

rn 

a 

tot.  mass  of  elastic  system  containing  the  m.s.a. 

<^a> 

t.a.  and  space-average  mean-square  velocity 

(of  elastic  system  containing  set  a)  due  to  all  m.s.a. 

Q  ro 
V 

t.a.  and  s.a.  mean  square  velocity  of  single  m.s.a. 

Sa 

s.a.  spectral  density  of  forces  acting  on  m.s.a. 

^a 

dissipation  loss  factor  for  set  a 

loss  factor  pertaining  to  power  flow  from  set  a  to  set  0 

% 

center  frequency  of  band  containing  all  modal  resonances 
of  both  sets 

set-to-set  power  flow  coefficient 

0o 

mode-to-mode  power  flow  coefficient  (averaged  over 
both  sets) 
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